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ABSTRACT 


This  report  documents  the  work  conducted  in  phase  I of  an  oixjoing  research 
program  of  which  the  c±>jective  is  to  develop  a head-neck  simulator  with 
omni-directional  biofidelity.  The  report  is  composed  of  two  volumes: 

Vol.  I : Analysis  volunteer  tests 
Vol.  II  : Mathematical  simulations 

In  Vol.  I a detailed  analysis  is  carried  out  of  a large  number  of  human 
volunteer  tests  conducted  by  the  Naval  Biodynamics  Laboratory  in  New 
Orleans.  Ihese  human  subjects  were  exposed  to  frontal,  lateral  and  oblique 
inpacts  with  an  impact  severity  up  to  15  g and  17  m/s.  This  analysis  re- 
sults for  each  impact  direction  in  a simple  analog  system  that  caqpletely 
specifies  the  observed  dynamical  behaviour. 

Vol.  II  documents  the  validation  of  the  proposed  analog  systems  and  com- 
pares the  behaviour  of  the  Part  572  and  Hybrid  III  head-neck  systems  with 
the  human  volunteer  behaviour.  It  follows  that  both  neck  designs  are  stif- 
fer  than  any  of  the  volunteers  in  the  impact  ranges  tested.  Preliminary 
results  are  presented  of  ccmputer  aided  design  activities  in  order  to 
inprove  existing  designs  in  view  of  these  new  findings. 
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CHAmR  1 
INTRODUCTION 

1. 1 Project  description  and  objectives 

The  National  Highway  Traffic  Safety  Administration  (NHTSA)  has  available 
large  quantities  of  test  data  on  the  head  and  neck  response  of  human  volun- 
teer subjects  in  frontal,  oblique  and  lateral  sled  tests.  The  Vehicle 
Research  and  Test  Center  (VRTC)  will  use  this  data  for  the  development  of 
an  inproved  head/neck  motion  simulator.  Where  possible,  this  human  test 
data  will  be  supplemented  by  data  of  human  cadavers  tested  at  exposure 
levels  beyond  those  generally  acceptable  for  human  volunteers. 

One  or  more  candidate  mechanical  systems  will  be  identified  to  reproduce 
observed  head/neck  motions.  Existing  models  will  form  the  basis  for  this 
analysis  but  the  study  will  not  necessarily  be  limited  to  such  models.  The 
performance  of  candidate  mechanical  neck  simulators  will  be  analyzed  and 
optimized  using  mathematical  model  simulations  resulting  finally  in  design 
specification  for  an  inproved  dummy  neck. 

This  report  presents  the  results  of  phase  I of  this  ongoing  research  pro- 
gram. The  report  is  conposed  of  two  volumes: 

Volume  I : Analysis  volunteer  tests 
Volume  II:  Mathematical  simulations. 

In  Vol.  I an  analysis  is  presented  of  a number  of  human  volunteer  tests  in 
three  different  impact  directions  conducted  at  the  Naval  Biodynamics  Labo- 
ratory (NBDL)*.  This  analysis  results  for  each  irtpact  direction  in  an 
analog  system  that  corpletely  specifies  the  observed  dynamical  behaviour. 

Vol.  II  documents  the  validation  of  the  prop<3sed  analog  systems  by  means  of 
mathanatical  simulations  and  presents  a number  of  preliminary  corputer 
aided  design  activities  in  order  to  improve  existing  dumnny  necks  in  view  of 
these  new  findings. 

* NBDL:  Naval  Biodynamics  Laboratory  in  New  Orleans;  Previously  Naval  Aero- 
space  Medical  Research  Laboratory  Detachment  (NAMRLD). 
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1 . 2 Summary  research  progran:  analyses  volunt-ooi  t.osts 

Realistic  simulation  of  the  neck  response  in  a dummy  is  of  vital  importance 
to  get  a humanlike  dynamical  behaviour  of  the  head.  Trajectories  of  the 
head  and  the  nature  of  the  head  contact  with  vehicle  interior  or  exterior 
are  critically  dependent  on  the  dummy's  neck  design.  In  literature  neck 
performance  criteria  are  limited  to  the  neck  response  in  forward  flexion 
and  extension.  Section  1.3  will  review  these  requirements  and  it  will  be 
shown  that  the  present  requirements  are  not  sufficient  conditions  to  ensure 
humanlike  response  in  frontal  impact  direction.  Moreover,  increasing  re- 
search activities  in  the  field  of  side  impact  and  pedestrian  protection 
clearly  show  the  need  for  additional  requirements  in  lateral  and  oblique 
direction. 

In  the  past  years  the  Naval  Biodynamics  Laboratory  (NBDL)  has  conducted  an 
extensive  research  program  to  determine  the  dynamic  head-neck  response  of 
volunteer  subjects  under  various  test  conditions  and  impact  directions.  In 
section  1.4  a brief  review  will  be  presented  of  relevant  literature  in  this 
respect.  Up  to  now  no  concrete  reccmmendations  for  performance  requirements 
have  been  based  on  these  NBDL  tests. 

The  analysis  of  the  human  volunteer  tests  in  this  study  is  based  on  de- 
tailed test  results  on  magnetic  tape,  obtained  frcm  NBDL.  Chapter  2 will 
summarize  the  NBDL  test  methodology,  test  conditions  and  the  instrumen- 
tation used  in  the  tests.  A description  will  be  given  of  the  type  of  test 
parameters  available  on  the  magnetic  tapes  and  of  the  selected  test  data 
used  in  our  study.  Coordinate  system  definitions  by  NBDL  anc]  the  human 
subject  antroponetry  will  be  described  in  chapter  3. 

A computer  program  was  developed  that  allows  visualisation  of  the  head  and 
Tl  displacements  on  the  basis  of  the  test  results  available  on  magnetic 
tape.  Chapter  4 will  present  results  for  the  Tl  displacements  and  chapter  5 
for  the  relative  head  motion.  Significant  variations  can  be  observed  in  the 
initial  neck  length  in  different  tests  with  the  same  subject  (see  chapter 
4).  It  will  be  shown  that  these  variations  are  mainly  due  to  errors  in  the 
specification  of  the  Tl  coordinate  system.  A method  will  be  proposed  to 
correct  for  these  errors. 
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From  the  analysis  of  the  relative  head  displacements  in  chapter  5 it  fol- 
lows that  the  motions  in  frontal,  lateral  as  well  as  oblique  impacts  can  be 
quite  realistically  approximated  by  means  of  a 2-pivot  linkage  with  2 
degrees  of  freedom  in  frontal  impacts  and  3 degrees  of  freedcm  in  lateral 
and  oblique  impacts. 

Chapter  6 deals  with  the  calculation  of  loads  applied  by  the  neck  to  the 
head  near  the  occipital  condyles  and  the  Tl  origin.  For  frontal  impacts, 
results  of  these  calculations  are  compared  with  the  response  corridors 
proposed  by  Mertz  et  al.  On  the  basis  of  these  load  calculations  dynamical 
elements  (joint  properties)  for  the  2-pivot  linkage  will  be  defined  (chap- 
ter 7).  Chapter  8 briefly  discusses  an  evaluation  of  sane  of  the  earlier 
tests  conducted  by  NBDL  and  a comparison  with  our  findings. 

The  major  findings  and  conclusions  of  this  research  program  as  far  as  the 
analysis  of  the  human  volunteer  tests  is  concerned,  are  discussed  and 
summarized  in  Chapter  9. 

1.3  Literature  review:  performance  requirements 

In  the  early  seventies  performance  requirements  for  dynamical  fiead-neck 
behaviour  in  frontal  impacts  were  proposed  by  Mertz  et  al.  (1,2),  based  on 
static  and  dynamical  tests  with  volunteers  and  human  cadavers.  Neck  re- 
sponse was  defined  by  the  relation  between  the  torque  about  the  occipital 
condyles  and  the  angular  position  of  the  head.  Resulting  performance  corri- 
dors specifically  for  the  loading  phase  are  shown  in  Fig.  1.  For  unloading 
motion  a ratio  was  defined  for  the  area  under  the  unloading  curve  with 
respect  to  the  area  under  the  loading  curve.  Separate  corridors  for  flexion 
and  extension  motions  were  developed.  Human  volunteer  dynamical  tests  in 
this  study  were  limited  to  one  subject,  partly  tested  with  an  additional 
weight  attached  to  the  head.  An  evaluation  of  several  neck  designs  with 
respect  to  these  requirements  is  incorporated  in  ref. (2). 


MOMENT  ABOUT  OCCIPITAL  CONDYLES 
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Flexion 


Fig.  i Head-neck  response  envelops  for  the  loading 
phase  according  to  Mertz  et.  al.  (2). 


Such  requirements,  however,  are  not  sufficient  conditions  to  ensure  human- 
like response  (2).  Both  Mertz  et  al.  (2)  and  Melvin  et  al.  (3)  discuss 
additional  displacement  requirements  such  as  a center  of  gravity  trajec- 
tory, but  concrete  recanmendations  in  this  respect  have  not  teen  formu- 
lated. 

A more  recent  paper  by  Muzzy  and  Lustick  (4)  illustrates  the  need  for  these 
type  of  requirements.  In  their  study  displacements  of  the  Hybrid  II  head- 
neck  system  were  compared  with  human  volunteers  in  10  g frontal  impacts. 
The  complete  dummy  was  tested  in  a restraint  system  similar  to  that  of  the 
human  volunteers.  Significant  differences  in  the  displacements  of  the  head 
center  of  gravity  and  the  head  angle  relative  to  the  torso  were  observed: 
the  Hybrid  II  neck  was  found  to  be  stiffer  than  the  volunteer. 

A preliminary  torque  rotation  response  corridor  for  lateral  flexion  was 
proposed  by  Patrick  and  Chou  (5).  This  corridor  which  is  shown  in  Fig.  2, 
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was  based  on  4 dynamic  tests  with  2 human  volunteers.  Torque-rotation 
curves  with  respect  to  the  occipital  condyles  for  the  individual  fests  are 
included  in  Fig.  2.  Considerable  differences  between  Ix^th  volunteers  can  !)e 
observed . 


SAK  15h  (K.II),  12.9  mph.  6.7g) 

SAE  1^5  (KID,  ID. A mpli, 


(AO, 40)  (50,401 


Fig.  2 A preliminary  torque  rotation  response  corridor  for 
lateral  flexion  according  to  Patrick  and  Chou  (5). 


Neck  perfonnance  requirements  used  in  a recent  evaluation  of  four  siuo 
impact  dummies  were  limited  to  the  maximum  head  lateral  flexion  as  observer; 
in  volunteer  and  cadaver  tests  neglecting  the  typical  three-dimens ion a J 
nature  of  lateral  head  neck  motions  in  a dynamical  environment  (6). 

In  contrast  to  frontal  and  lateral  impacts,  for  oblique  motion  no  response 
requirements  are  reported  in  literature. 
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1.4  Literature  review:  human  volunteer  tests  conducted  by  NBDL 

A detailed  description  of  the  test  methodology  in  the  NBDL  tests  is  pro- 
vided in  references  (7-12)  and  will  be  summarized  in  section  2.1  (test 
method).  The  first  experiments  in  the  NBDL  test  program  were  conducted  on 
the  WHAM  II  accelerator  at  Wayne  State  University  in  the  period  between 
1967  and  1969.  In  total  236  experiments  with  15  subjects  in  frontal  flex- 
ion were  conducted.  Head  and  T1  accelerations  and  velocity  time  histories 
of  18  representative  runs  involving  6 subjects  are  presented  in  references 
(13-15).  More  detailed  results  including  head  displacement- time  histories 
with  respect  to  Tl,  can  be  found  in  reference  (7).  As  a part  of  our  study 
an  analysis  has  been  made  of  displacement  data  provided  in  this  reference 
(see  chapter  8). 

By  Ewing  and  Thcanas  (16)  torque-angular  displacement  characteristics  for 
seme  of  these  tests  were  detennined  and  compared  with  data  of  Mertz  et  al. 
(1,2).  EXie  to  different  instrumentation  and  measuring  techniques  several 
differences  were  noted,  in  particular  in  the  initial  phase  of  the  motion 
where  the  sign  of  the  occipital  torque  was  the  opposite  of  the  Mertz  data. 
For  larger  head  excursions  torque  data  were  found  to  be  conparable. 

A subsequent  series  of  100  frontal  iirpact  runs  with  13  volunteers  was  con- 
ducted in  the  early  seventies  on  the  HYGE  accelerator  of  the  laboratory  of 
NBDL.  Newly  developed  instrumentation  and  photographic  mount  were  used.  The 
main  objective  of.  this  study  was  the  analysis  of  the  effect  of  the  initial 
head  position  on  head  linear  and  angular  accelerations  and  angular  velo- 
cities (17),  A considerable  influence  of  the  head-neck  initial  position  was 
reported.  In  an  additional  test  program  of  75  tests,  the  effect  of  sled 
pulse  duration,  rate  of  onset  and  peak  acceleration  on  head  accelerations 
and  head  velocities  was  evaluated  (18).  The  experimental  results  were  found 
to  be  affected  considerably  by  these  changes  in  experimental  conditions. 

The  pelvic  response  in  these  type  of  impacts  was  measured  in  a test  program 
with  one  subject.  Significant  differences  were  observed  between  Tl  and 
pelvis  acceleration-time  histories  due  to  relatively  larger  motion  of  Tl 
permitted  by  the  restraint  system  and  due  to  chest  deformation  (19). 
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More  recently  a large  number  of  lateral  and  oblique  tests  were  conducted  by 
NBDL.  A detailed  description  of  test  method  and  test  results  for  the  first 
serie  of  lateral  tests  (i.e.  34  experiments  with  6 subjects)  is  provided  iti 
ref.  (10).  Based  on  additional  test  data  obtained  from  NBDL,  in  our  study 
an  analysis  was  made  of  some  of  the  most  severe  tests  in  this  test  program 
(see  chapter  8). 

Ref.  (20)  presents  for  a subsequent  test  series  with  5 subjects,  head 
accelerations  and  velocities  as  function  of  time.  Various  sled  acceleration 
profiles  were  used  with  peak  values  up  to  11  g. 

In  ref.  (21)  the  effect  of  initial  head-neck  position  on  head  accelera- 
tions, velocities  and  displacements  for  100  lateral  tests  with  6 human 
subjects,  was  evaluated.  Peak  sled  accelerations  in  these  experiments 
ranged  fron  2-7  g.  As  with  for  the  frontal  impacts  a significant  effect  of 
the  initial  head-neck  position  was  observed.  One  of  the  subjects  incorpo- 
rated in  this  test  program  i.e.  subject  H00083,  is  also  incorporated  in  the 
database  to  be  analysed  in  the  present  study. 


j 

I 
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CflAlTtR  2 

DATABASE  DESCRIPl’ION  (SUBJECTS  H00083  ANO  H0009  3) 

2. 1 Test  method 

In  the  tests  conducted  at  NBDL  the  subjects  are  seated  on  a Bendix  HYGE 
pneumatically  driven  0.3048  diameter  accelerator.  The  sled  mass  is  1669  kg 
and  the  acceleration  stroke  is  limited  to  1.52  m (12).  In  the  frontal 

impact  tests  the  subjects  are  in  an  upright  position  restrained  by  shoulder 
straps,  a lap  belt  and  an  inverted  V-pelvic  strap  tied  to  the  lap  belt. 
Upper  arm  and  wrist  restraints  are  used  to  prevent  flailing  (12).  In  addi- 
tion a loose  safety  belt  around  the  chest  is  employed.  The  same  restraint 
system  is  used  in  lateral  and  oblique  tests.  In  addition  a lightly  padded 
wooden  board  is  placed  against  the  right  shoulder  of  the  subject  to  limit 
the  upper  torso  motion. 

In  frontal  impacts  the  thrust  vector  is  ncminally  directed  from  chest  to 
back  and  in  lateral  inpacts  fran  the  right  to  the  left  shoulder  (Fig.  7). 
In  oblique  tests  the  thrust  vector  is  directed  fron  right-front  to  left- 
back  (angle  between  thrust  vector  and  left-right  axis  is  45  degrees). 

2.2  Instrumentation 

The  resulting  three-dimensional  motions  of  the  head  and  first  thoracic 
vertebral  body  (Tl)  are  monitored  by  anatcxnically  mounted  clusters  of 
acceleroneters  and  photographic  targets.  Both  the  head  and  Tl  are  equipped 
with  six  piezo-resistive  acceleroneters  mounted  on  a T-shaped  plate  at  the 
mouth  and  at  the  spinous  process  of  Tl  (Fig.  3).  The  configuration  of  tlie 
accelerometers  and  the  error  propogation  associated  with  this  methoci  for 
determining  kinematic  quantities  has  been  described  by  Becker  and  Willems 
(9).  More  recently  the  application  of  a nine-accele rone ter  configurations 
at  NBDL  has  also  been  reported  (22). 
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Fig.  3 Six-accele rone ter  system  3-2-1  configuration  (12). 


In  addition  to  the  accelerometers,  a two  axis  rate  gyroscope  is  mounted  on 
both  T plates.  The  resulting  data  are  used  as  an  independent  measurement 
of  two  components  of  the  angular  velocity  of  the  head  and  neck  (12).  Test 
data,  obtained  within  the  present  study  fron  NBDL,  did  not  include  this 
type  of  gyroscopic  based  data. 

Groups  of  black  and  white  phototargets  are  connected  to  the  T plates  car- 
ried by  head  and  Tl  (Fig.  4).  Cinephotographic  coverage  of  these  targets  is 
provided  by  sled  mounted  camera's  (pin-registered  16  mm,  500  frames  per 
second),  situated  orthogonally  to  each  other  at  approximately  one  meter 
from  the  subject  (8,12).  This  method  permits  three-dimensional  displace- 
ments of  head  and  Tl  to  be  determined  relative  to  the  sled. 

In  order  to  compare  subjects  at  similar  points  in  the  anatcmy,  it  is  neces- 
sary to  define  head  and  Tl  anatomical  coordinate  system  (see  chapter  3). 
These  antropometric  coordinate  systems  are  related  to  instrumentation 
coordinate  systems  by  means  of  three-dimensional  X-ray  technigues  (11). 
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Fig.  4 Instrumented  human  volunteer  (9). 

2.3  Processed  data  obtained  from  NBDL 

Data  processing  by  NBDL  includes  conversion  of  photo  and  sensor  data  into 
kinematics  of  the  head  and  T1  anatonical  coordinate  systens.  Ihe  solution 
for  the  kinematic  variables  derived  frcm  acceleraneter  data  requires  ini- 
tialisation of  location  and  orientation  frcm  photo  derived  data  and  inte- 
gration of  the  non-linear  differential  equations  defined  by  expressions  for 
the  acceleraneter  measurements  (9).  The  solution  for  displacements  and 
angular  rotations  is  obtained  by  a least  square  algorithm,  which  makes  use 
of  measurements  in  the  film  plane  frcm  two  different  cameras  (8).  Succes- 
sive differentiation  is  used  to  determine  velocity  conponents  fron  the 
photo  derived  displacements. 

In  this  way  for  each  test  91  kinematic  variables  as  function  of  time  are 
available  on  the  tapes  obtained  frcm  NBDL.  Table  1 provides  a summary  of 
this  data.  Conponents  of  rotations,  translations,  linear  velocities  and 
linear  accelerations  are  all  expressed  with  respect  to  the  laboratory  or 
sled  coordinate  system  while  conponents  of  angular  velocities  and  accele- 
rations are  given  with  respect  to  anatomical  coordinate  systems. 
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Table  1 - Test  results  for  the  head,  T1  and  sled  motions 

available  on  the  NBDL  data  tapes 

- Translations  and  linear  velocities  of  the  head  and  T1 
anatcraical  origins  with  respect  to  the  sled,  derived 
fran  photographic  as  well  as  acceleroneter  data  (24)*. 

- Translations,  linear  velocities  and  linear  accelerations 
of  head  and  T1  anatomical  origins  with  respect  to  the 
laboratory,  derived  from  accelerometer  data  (18). 

- Angular  velocities  of  the  head  and  Tl  about  the  anatomical 
coordinate  axes  derived  from  photographic  as  well  as  acce- 
lercsmeter  data  (12). 

- Angular  accelerations  of  the  head  and  Tl  about  the  anato- 
mical coordinate  axes,  derived  from  accelerometer  data  (6), 

- Rotations  of  the  head  and  Tl  expressed  in  euler  angles  as 
well  as  quaternians  with  respect  to  the  laboratory,  derived 
from  photographic  as  well  as  accelerometer  data  (28). 

- Acceleration,  velocity  and  displacement  of  the  sled  with 
respect  to  the  x-axis  of  the  laboratory  (3). 


* Number  of  variables  in  this  category. 


Some  of  the  kinematic  variables  are  derived  independently  from  accelero- 
meter as  well  as  photographic  measurements  (i.e.  linear  translations  and 
linear  velocities  relative  to  the  sled,  rotations  and  angular  velocities). 
A coTparison  between  these  data  showed  discrepancies  in  a nuimber  of  cases 
particularly  in  some  of  the  tests  with  large  three-dimensional  motions 
(head  motion  in  lateral  and  oblique  impacts).  Such  differences  also  have 
been  reported  by  NBDL  (23).  Optimal  combination  of  the  results  from  both 
sets  of  measurements  into  one  consistent  set  of  kinematic  variables  is  one 
of  the  objects  of  the  NBDL  research  program  (23). 


Our  further  analysis  will  be  based  on  a subset  out  of  the  91  NBDL  vari- 
ables. Table  2 summarizes  these  variables  together  with  a code  defined  by 
NBDL  to  identify  these  variables  on  the  tapes.  A coiplete  description  of 
all  91  variables  is  provided  in  ref.  (10).  The  displacements  data  (trans- 
lations and  rotations)  in  this  subset  are  derived  from  photographic  data 
while  the  accelerations  are  derived  from  accelerometer  data.  Head  angular 
velocities  which  are  required  for  the  load  calculations  described  in 
chapter  6,  are  based  on  accelercroeter  data. 
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Table  2 Selected  set  of  kinematic  variables  on  NBDf,  data  ta^x?s 


NBDL 

Definition 

Description 

AAXXOS 

AAYXOS 

AAZXOS 

Linear  accelerations  of  the  head  anataxiical  origin  along  the 
x-axis  (AAXXOS),  the  y-axis  (AAYXOS)  and  the  z-axis  (AAZXOS) 
respectively  of  the  laboratory  coordinate  system  with  respect 
to  the  fixed  laboratory  coordinate  system  as  derived  frcxn 
accelerometer  data. 

ANXXOS 

ANYXOS 

ANZXOS 

Linear  accelerations  of  the  T1  anatomical  origin  defined 
corresponding  to  the  linear  head  accelerations  AAXXOS,  AAYXOS 
and  AAZXOS. 

QHAOXS 

OHBOXS 

OHCOXS 

Angular  accelerations  of  the  head  about  the  x-axis  (QHAOXS), 
y-axis  (QHBOXS)  and  z-axis  (QHCOXS)  respectively  of  the  heaa 
anatomical  coordinate  system  as  derived  fron  acceleraneter  data 

ONAOXS 

QNBOXS 

QNCOXS 

Angular  accelerations  of  Tl  about  the  x-axis  (QNAOXS),  y-axis 
(QNBOXS)  and  z-axis  (QNCOXS)  respectively  of  the  Tl  anatomical 
coordinate  system  as  derived  frcm  acceleraneter  data. 

RhiAOXS 

RH130XS 

RHCOXS 

Angular  velocities  of  the  head  about  the  x-axis  (RHAOXS),  the 
y-axis  (RHBOXS)  and  the  z-axis  (RHCOXS)  respectively  of  the 
head  anatomical  coordinate  system  as  derived  fron  accelero- 
meter data. 

ACXXOS 

VCXXOS 

Linear  acceleration  (ACXXOS)  and  linear  velocity  (VCXXOS)  of 
the  sled  along  the  x-axis  of  the  laboratory  coordinate  system 
with  respect  to  the  fixed  laboratory  coordinate  systems  as 
derived  from  sled  mounted  acceleraneter  data. 

DAXSOP 

Qf^YSOP 

DAZSOP 

Linear  displacements  of  the  head  anatomical  origin  along  the 
x-axis  (DAXSOP),  the  y-axis  (DAYSOP)  and  the  z-axis  (DAZSOP) 
respectively  of  the  sled  coordinate  system  with  respect  to  the 
sled  coordinate  system  as  derived  from  photo  target  data. 

DNXSOP 

DNYSOP 

DNZSOP 

Linear  displacements  of  the  Tl  anatonical  origin  defined 
corresponding  to  the  linear  head  displacements  DAXSOP,  DAYSOP 
and  DAZSOP. 

4H001P 

4H002P 

4H003P 

4H004P 

Qaarternians  - The  four  variables  which  define  the  angular 
orientation  of  the  head  anatonical  coordinate  system  relative 
to  the  laboratory  coordinate  system  as  derived  fron  photo 
target  data  (see  eg.  (1)). 

4N001P 

4N002P 

4N003P 

4N004P 

Quarternians  - The  four  variables  which  define  the  angular 
orientation  of  the  Tl  anatomical  coordinate  system  relative  to 
the  laboratory  coordinate  system  as  derived  fron  photo  target 
data. 

PHC03P 

Angle  of  rotation  (euler  angle)  of  the  head  about  the  z-axis  of 
the  head  anatonical  coordinate  system  as  derived  from  photo 
target.  (Head  anatonical  coordinate  systan  is  initially 
aligned  with  the  laboratory  coordinate  system). 
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Qjarternians  will  be  used  to  define  the  angular  orientation  head  and  Tl. 
Tlie  transformation  matrix  to  transfocm  a vector  from  ttie  anatomical  co»,:)r- 
dinate  system  to  the  laboratory  coordinate  system  in  terms  of  quarternians 
is  as  follows  (10): 


In  addition  to  quarternians  euler  angles  are  also  provided  on  tape.  These 
angles  are  defined  by  NBDL  as  follows:  The  angular  orientation  of  a body  is 
thought  to  be  the  result  of  three  successive  rotations.  Before  the  first 
rotation  the  coordinate  system  of  the  body  (head  or  Tl)  coincides  with  the 
laboratory  (sled)  coordinate  system.  The  first  rotation  is  carried  out 
about  the  x-axis.  The  second  rotation  is  performed  about  the  carried  y-axis 
of  the  body  and  the  third  rotation  about  the  carried  z-axis.  The  third 
euler  angle  for  the  head  (denoted  by  PHC03P  in  Table  2)  will  be  used  here 
for  evaluation  of  the  head  twist  in  lateral  and  oblique  irrpacts. 

Accelerometer  derived  data  are  provided  on  the  tapes  with  a time  interval 
of  0.0005  s.  E^otographic  derived  variables  are  related  to  the  time-base  by 
a separate  time-channel  (time  interval  ^ 0.002  s). 

2. 4 Testconditions  of  30  selected  tests 

Datatapes  containing  the  results  of  85  human  volunteer  tests  with  subject 
H00083  and  H00093  were  obtained  from  NBDL.  The  number  of  frontal,  lateral 
and  oblique  tests  in  this  database  is  14,  30  and  41  respectively.  Out  of 
these  tests  a subset  of  30  of  the  most  severe  experiments  is  selected  to  be 
used  for  the  analysis  presented  in  this  report.  Table  3 summarizes  these 
tests  and  the  most  important  test  characteristics.  The  selection  of  this 
subset  was  based  on  the  following  criteria: 


2 2 2 2 
Q4  + Ql  - Q2  - qg 


A = 2(q^q^  + q^q2) 

2(qiq3  " P4q2) 


where:  q^  = 4H001P  (4N001P) 


= 4H002P  (4N002P) 

q^  = 4H003P  (4N003P) 

Q = 4H004P  (4N004P) 
^4 
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- no  large  data  gaps  of  data  shifts  iii  the  photcxjrapliic  dot  ived  varia[)l>.‘.s 

- emission  of  test  with  a low  impact  level  (i.e.  low  sle<J  velocities  and/or 
low  sled  accelerations) 

- no  large  forward  or  lateral  bending  of  the  head  and  neck  in  the  initial 
position. 


Table  3 - Test  characteristics  of  30  selected  runs  with  2 human 


volunteers  (H00083  and  H00093). 


Subject 

Test 

number 

Sled  velocity 
change  (m/s) 

Peak  sled2 
acc.  (m/s  ) 

Rate^of  onset* 
(m/s  ) 

FRONTAL 

H00083 

LX3530 

7.5 

49 

982 

H00083 

LX3536 

10.5 

73 

1613 

H00083 

LX3544 

12.4 

92 

2205 

H00093 

LX3531 

7.6 

51 

1050 

H00093 

LX3537 

10.1 

71 

1585 

H00093 

LX3548 

12.2 

90 

2142 

H00093 

LX3550 

13.3 

98 

2498 

H00093 

LX3558 

14.3 

110 

2903 

H00093 

LX3578 

16.0 

132 

3898 

H00093 

LX3583 

16.7 

142 

4269 

H00093 

LX3616 

17.2 

150 

4873 

LATERAL 

H00083 

LX1831 

6.5 

49 

3229 

H00083 

LX2013 

6.3 

49 

1079 

H00083 

LX2124 

3.1 

89 

9435 

H00083 

LX2302 

6.4 

69 

1497 

H00093 

LX2151 

3.3 

100 

11075 

H00093 

LX2326 

6.5 

60 

5764 

H00093 

LX2355 

6.5 

50 

4961 

OBLIQUE 

H00083 

LX2786 

7.5 

50 

7061 

H00083 

LX2872 

9.0 

61 

8642 

H00083 

LX3093 

10.3 

71 

1595 

H00083 

LX3102 

9.0 

60 

1337 

H00083 

LX3133 

8.8 

79 

1839 

H00083 

LX3153 

8.8 

80 

8461 

H00093 

LX2784 

7.5 

50 

6802 

H00093 

LX2843 

3.0 

89 

1302 

H00093 

LX3122 

8.9 

79 

1857 

H00093 

LX3148 

8.9 

97 

2573 

H00093 

LX3153 

8.8 

78 

8669 

H00093 

LX3417 

8.8 

59 

3404 

* Slope  of 

the  best 

least  square  line 

fit  of  the  rising 

portion  of  the  sled 

acceleration  profile  between  20%  and  50%  of  the  peak  sled  acceleration. 


1 


16 


The  resulting  subset  contains  11  frontal,  7 lateral  and  12  oblique  iiTpacts. 
Sled  acceleration  time  histories  (ACXXOS)  for  the  tests  in  this  database 
are  incorporated  in  Annex  A.  Fig.  5 presents  a comparison  of  the  sled 
accelerations  in  the  most  severe  frontal,  lateral  and  oblique  test  (i.e. 
LX3616,  LX2302  and  LX3148,  respectively).  The  oblique  test  appears  to  be 
more  severe  than  the  lateral  one,  while  the  frontal  test  is  more  severe 
than  the  oblique  one. 


Fig.  5 Comparison  of  sled  accelerations  in  the  most 
severe  frontal,  lateral  and  oblique  test. 
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CHAPTER  3 

COORDINATE  SYSTEMS  AND  HUMAN  SUBJ^XT  ANrHROPOMEl’KY 
3.1  Introduction 


Definitions  taken  for  the  coordinate  systene  are  identical  to  those  used  by 
NBDL  (7,10,17).  Section  3.2  will  discuss  the  head  and  Tl  anatcmical  coor- 
dinate systonns  and  section  3.3  the  laboratory  and  sled  coordinate  systems. 
In  addition  instrumentation  and  camera  coordinate  systems  are  also  defined 
by  NBDL  but  these  systems  are  not  necessary  to  understand  the  results 
presented  in  this  report. 

In  order  to  calculate  neck  loads,  estimations  have  to  be  made  for  the  head 
mass  distribution  and  the  location  of  the  occipital  condyles.  Partly,  these 
estimations  will  be  based  on  available  anthropometric  data  for  each  subject 
(section  3.4).  Section  3.5  will  present  estimations  for  the  head  mass  and 
moments  of  inertia  while  section  3.6  discusses  the  location  of  the  occi- 
pital condyles  and  the  center  of  gravity. 

3.2  Head  and  Tl  anatcmical  coordinate  systems 

The  head  and  Tl  anatcmical  coordinate  systems  defined  by  NBDL  are  orthogo- 
nal and  right-handed  (Fig.  6).  For  the  head,  the  anatomical  origin  is  at 
the  midpoint  of  the  line  connecting  the  superior  edges  of  the  auditor^^ 
meati.  The  positive  x-axis  is  the  axis  from  the  origin  through  the  midpoint 
of  the  line  connecting  the  infraorbital  notches.  The  positive  z-axis  is 
perpendicular  to  the  x-axis  in  superior  direction  in  such  a way  that  the 
xz-plane  is  the  midsagittal  plane.  Ihe  positive  y-axis  is  perpendicular  to 
the  xz-plane  toward  the  left  ear.  In  this  way  the  xy-plane  is  approximately 
the  Frankfort  plane. 
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Fig.  6 Location  of  head  and  Tl  anatomical  coordinate  systems 
(o.c.  = occipital  condyles,  c.g.  = center  of  gravity). 

The  origin  of  the  spine  (Tl)  anatomical  coordinate  system  is  at  the  ante- 
rior superior  comer  of  the  first  thoracic  vertebral  body  (Tl).  The  posi- 
tive x-axis  is  directed  frcm  the  most  posterior  point  of  the  spinous  pro- 
cess of  Tl  to  the  origin,  the  positive  z-axis  is  perpendicular  to  the 
x-axis  in  superior  direction  and  the  positive  y-axis  is  toward  the  left, 
perpendicular  to  the  xz-plane. 

Ihree-dimensional  X-ray  techniques  were  used  to  specify  these  anatomical 
coordinate  systems  relative  to  the  anatcmical  landmarks,  where  for  the 
identification  of  the  landmarks  in  the  head  small  lead  markers  were  located 
near  the  auditory  meati  and  the  infraorbital  notches.  The  location  of  the 
instrumentation  packages  with  respect  to  the  anatcmical  coordinate  systems 
were  determined  using  the  same  X-ray  measurements. 

3. 3 Laboratory  and  sled  coordinate  systems 

The  basic  reference  frame  in  the  NBDL  tests  is  the  laboratory  coordinate 
system,  which  is  a right-handed  orthogonal  coordinate  system  fixed  in  the 
laboratory  (Fig.  7).  The  positive  x-axis  of  this  system  is  defined  parallel 
but  in  opposite  direction  to  the  thrust  vector  of  the  accelerator.  The 
positive  z-axis  is  opposite  in  direction  to  gravity.  To  the  sled,  a sled 
coordinate  system  is  connected  with  the  same  orientation  as  the  laboratory 
coordinate  system. 
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Fig.  7 Itough  orientation  of  head  and  T1  anatonical  coordinate 
systems  relative  to  laboratory  and  sled  coordinate 
system  in  case  of  frontal,  oblique  and  lateral  impacts. 

Ihe  initial  naninal  orientation  (i.e.  before  moving  of  the  sled)  of  the 
head  and  T1  anatomical  coordinate  systems  relative  to  the  laboratory  and 
the  sled  coordinate  system  for  the  3 impact  directions,  is  roughly  indi- 
cated in  Fig.  7.  It  follows  that  the  naninal  value  of  the  third  euler  angle 
(i.e.  rotation  about  head  or  T1  z-axis  as  specified  in  section  2.3)  is  0 
degrees  in  frontal,  45  degrees  in  oblique  and  close  to  90  degrees  in  case 
of  lateral  impacts. 

Later  on  in  chapter  3 and  4 information  on  the  actual  initial  orientation 
of  Tl  and  head  anatomical  coordinate  systems  will  be  provided.  It  will  be 
shown  that  the  initial  oientation  of  the  Tl  anatonical  coordinate  system 
varies  considerably  in  different  tests  for  the  same  subject,  which  is  ex- 
plained by  errors  in  the  specification  of  this  coordinate  system  (chapter 
4). 
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3.4  Human  subject  anthropcynetric  measurements 


At  MBDL  detailed  anthropometric  measurements  are  conducted  for  each  subject 
as  part  of  the  test  protocol.  Table  4 summarizes  the  most  significant  data 
in  this  respect.  Definitions  for  these  ant ropome trie  variables  can  be  found 
in  ref . ( 7 ) . 


Table  4 - Selected  Human  Subject  Anthropometry 


Subject 

Number 

Standing 

Height 

(cm) 

Weight 

(kq) 

Sitting 

Height 

(cm) 

Head 

Circumference 

(cm) 

Head 

Breadth 

(cm) 

Head 

Length 

(cm) 

H00083 

180.0 

75.7 

92.4 

54.7 

15.5 

18.3 

H00093 

173.0 

72.5 

93.5 

52.1 

14.4 

17.8 

3. 5 Head  mass  and  moments  of  inertia 

Estimations  for  the  subjects  head  mass  and  the  moments  of  inertia  about  the 
principal  axes  will  be  based  on  regression  equations  proposed  by  McConville 
et  al.  (24).  Separate  equations  were  given  for  estimations  based  on  stature 
and  weight  and  based  on  head  circumference,  head  breadth  and  head  length. 
Since  estimations  using  head  anthropemetry  data  are  more  reliable,  these 
equations  will  be  applied  here.  The  following  set  of  regression  equations 
was  used  in  our  study: 


M = d(C^ 

X HC  - C 

2 

X 

HL  - C3) 

(2) 

I 

C 

X HC  + C 

X 

HB  - C 

(3) 

XX 

1 

2 

3 

I 

C 

X HC  - C 

(4) 

yy 

1 

3 

I 

C 

X HC  - C 

X 

HL  - C^ 

(5) 

zz 

1 

2 

3 

where ; 

M 

= the  head  mass  (kg) 

I 

= moment 

of 

inertia  about  principal  x-axis 

2 

(kgm  ) 

XX 

2, 

I 

= moment 

of 

inertia  about  principal  y-axis 

(kgm  ) 

yy 

. 2 

I 

= moment 

of 

inertia  about  principal  z-axis 

(kgm  ) 

zz 

3 

d = density  (kg/dm  ) 
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HC  = head  circumference  (m) 

HB  = head  breadth  (m) 

HL  = head  length  (m) 

Ihe  values  for  the  constants  and  are  summarized  in  Table  5. 


Table  5 

- Coefficients 

in  regression  eqs. 

(2-5). 

eq.(2)  : 

= 21.618 

= 12.184 

2 

= 5.5936 

eq. (3)  : 

C = 0.11639 

1 

C = 0.17924 

1 

C = 0.11857 

C = 0.10605 

2 

= 0.0625049 

eq . ( 4 ) : 
eq.(5)  ; 

C = 0.07320 

C = 0.0794181 
3 

C = 0.0382935 

1 

2 

3 

The  density  will  be  taken  1 here  based  on  recommendations  made  by 
McConville  et  al.  (24).  This  value  is  slightly  smaller  than  values  obtained 
fran  cadaver  segment  measurements,  to  account  for  reported  overestimations 
by  these  regression  equations  (24). 

Table  6 summarizes  the  masses  and  nonents  of  inertia  for  each  subject 
resulting  from  these  regression  equations.  These  values  have  to  be  cor- 
rected for  the  effect  of  the  head  instrumentation.  Based  on  data  presented 

by  Ewing  and  Thomas  (16)  the  correction  for  the  head  mass  is  + 0.53  kg  and 

2 

for  the  mcment  of  inertia  about  the  principal  y-axis:  + 0.0075  kgm  . A 

similar  correction  is  made  for  the  moment  of  inertia  about  the  principal 

x-axis,  while  the  correction  for  the  moment  of  inertia  about  the  principal 

2 

z-axis  is  estimated  to  be  very  small:  + 0.0005  kgm  . The  resulting  masses 
and  mcments  of  inertia  are  incorporated  in  Table  6. 
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Table  6 - Estimations  for  Masses  and  Manents  of  Inertia 


MI  I I 

XX2  yy2  2^2 

Subject  (kg)  (kgm  ) (kgm  ) (kgm  ) 


Without  Instrumentation 


H00083 

4.00 

0.0176 

0.0186 

0.0132 

H00093 

3.50 

0.0134 

0.0140 

0.0105 

Corrected 

for  Instrumentation 

H00003 

4.53 

0.0251 

0.0261 

0.0137 

H00093 

4.03 

0.0209 

0.0215 

0.011 

The  location  of  the  principal  inertia  axes  will  also  be  based  on  calcula- 
tions of  McConville  et  al.  (24):  the  principal  y-axis  is  taken  parallel  to 

the  anatonical  y-axis  and  the  angle  between  principal  z-axis  and  head 

o , 

anatomical  z-axis  is  taken:  -36  (principal  z-axis  is  rotated  backwards). 
These  values  are  close  to  data  reported  by  Beier  (25).  The  orientation  of 
the  principal  axis  is  assumed  here  not  to  be  affected  by  the  head  instru- 
mentation packages. 

3. 6 Location  of  center  of  gravity  and  occipital  condyles 

No  regression  equations  are  available  to  estimate  the  location  of  the 
center  of  gravity  and  the  location  of  the  occipital  condyles  as  a function 
of  anthropometric  data.  The  location  of  the  center  of  gravity  will  be  based 
here  on  average  values  reported  by  Beier  (25)  resulting  from  measurements 
on  21  fresh  cadaver  heads  (Table  7).  These  values  have  been  corrected  for 
the  head  instrumentation  based  on  data  by  Ewing  and  Thomas  (16).  The  re- 
sulting c.g.  location  is  incorporated  in  table  7. 

'The  location  of  the  occipital  condyles  will  be  based  on  average  values 
reported  by  Ewing  and  Thomas  (16)  resulting  frcm  measurements  on  12  human 
subjects  (Table  7).  The  occipital  condylar  point  is  defined  by  Ewing  and 
Thomas  as  the  point  most  representative  of  the  center  of  the  condyles  as 
seen  on  lateral  X-ray  and  projected  onto  the  midsagittal  plane  of  the  head. 
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The  approximate  location  of  the  center  of  gravity  (c.g.)  and  the  occipital 
condylar  point  (o.c. ) is  indicated  in  Fig.  6. 


Table  7 - Location  of  center  of  gravity  and  occipital  condyles 
with  respect  to  head  anatonical  coordinate  system 


X 

y 

z 

Center  of  gravity 

(cm) 

(cm) 

(cm) 

Beier  (25) 

0.83 

0 

3.12 

Shift  due  to  instrumentation (16) 

+0.35 

0 

-0.20 

Resulting  c.g.  location 

1.2 

0 

2.9 

Occipital  condylar  point 

-1.1 

0 

-2.6 
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CHAPTER  4 

ANALYSIS  OF  Tl  MOTION  AND  CORRECTION  OF  Tl  COORDINATE  SYSTET^ 

4.1  Introduction 

In  order  to  visualize  volunteer  test  results,  a coirputer  program  was  deve- 
loped which  generates  a graphic  representation  of  the  three-dimensional 
head  and  spine  displacements,  based  on  photographic  derived  data  on  the 
NBDL  tapes.  This  program  generates  two  types  of  plots:  The  first  plot 
contains  displacements  in  a vertical  plane  with  the  same  orientation  as  the 
laboratory  (x,z)-plane.  Since  this  plane  is  parallel  to  the  thrust  vector 
of  the  accelerator,  it  will  be  called  plane  of  impact.  The  second  plot 
contains  displacements  in  a vertical  plane  perpendicular  to  the  plane  of 
inpact  and  with  the  same  orientation  as  the  laboratory  (y,z)-plane. 

For  the  vertebra  Tl,  the  trajectory  of  the  Tl  anatomical  origin  relative  to 
the  sled  is  generated,  together  with  the  Tl  anatcmical  z-axis  at  specified 
time  frames.  Section  4.2  will  present  resulting  plots  (in  the  impact  plane) 
for  the  30  tests  in  the  selected  subset. 

Section  4.3  evaluates  the  initial  neck  length  which  is  defined  as  the 
distance  between  head  and  Tl  anatomical  origin.  This  quantity  shows  large 
variations  in  different  tests  with  the  same  subject  which  will  be  explained 
by  errors  in  the  specification  of  the  Tl  coordinate  system.  In  section  4.4 
a corrected  Tl  coordinate  system  will  be  proposed. 

4.2  Tl  motion  in  the  plane  of  impact  and  Tl  initial  orientation 

Tl  moticxis  in  the  plane  of  inpact  are  presented  in  Fig.  8.  For  all  three 
impact  directions,  trajectories  of  the  Tl  anatonical  origin  relative  to  the 
sled  and  the  corresponding  position  of  the  Tl  anatomical  z-axis  at  various 
time  frames,  are  shown.  The  relative  Tl  origin  displacements  in  horizontal 
direction  appear  to  vary  frcrn  0.04  m for  the  less  severe  up  to  0.075  m for 
the  more  severe  impacts. 
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Fig.  8 indicates  small  vertical  translations  and  small  rotations  of  the  Tl 
vertebra  in  response  to  the  iitpact.  Tliese  displaceir^nts  are  considered 
small  enough  to  justify  neglecting  them  further  in  our  analysis.  The  same 
holds  for  the  out  of  impact  plane  displacements  which  also  are  found  to  be 
very  small.  Consequently  the  only  significant  displacement  of  Tl  is  a 
horizontal  translation  in  the  inpact  direction.  Corresponding  Tl  accelera- 
tions as  function  of  time  together  with  the  sled  deceleration  are  presented 
in  Fig.  9 for  three  inpact  directions.  For  each  inpact  direction  results  of 
one  typical  severe  test  are  shown.  Annex  A presents  this  type  of  data  for 
all  30  selected  tests. 


— time  I ms) 


0 40  80  120  160  200  2U0  280 


Fig.  9 Typical  exarrples  of  horizontal  Tl  and  sled  acceleration 
as  function  of  time  in  three  different  impact  directions. 
Tl; = sled. 

In  frontal  inpacts  this  Tl  acceleration  curve,  except  for  seme  small  os- 
cillations and  a phase  shift,  appears  to  be  close  to  the  sled  acceleration 
curve.  In  lateral  and  oblique  inpacts,  however,  a much  larger  deviation 
between  both  curves  can  be  observed,  indicating  the  different  nature  of  the 
interaction  between  the  restraint  system  and  torso  or  shoulder  in  these 
type  of  tests.  A cenparison  of  the  initial  orientation  of  the  Tl  anatanical 
coordinate  system  in  different  tests  for  the  same  subject  shows  in  a number 
of  cases  considerable  variations.  This  is  clearly  illustrated  for  instance 
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for  tests  LX3558  and  LX3616  with  subject  H00093  in  frontal  flexion  (see 
Fig.  8).  Such  deviations  which  cannot  be  fully  explained  tran  variations  in 
the  initial  torso  orientation,  are  roost  likely  caused  by  errors  in  the 
specification  of  the  Tl  coordinate  system.  Section  4.4  will  present  a 
corrected  Tl  coordinate  system  of  which  the  initial  orientation  is  identi- 
cal to  the  orientation  of  the  sled  (laboratory)  coordinate  system. 

4. 3 The  initial  neck  length 

The  initial  neck  length  will  be  defined  here  as  the  distance  between  Tl  and 
head  anatomical  origin  at  time  = 0 (i.e.  before  moving  of  the  sled).  Fig. 
10  presents  a ccstrparison  of  this  parameter  for  all  tests  obtained  fron  NBDL 
(i.e.  85  experiments).  For  the  same  subject,  variations  up  to  0.07  m in 
this  quantity  can  be  observed. 

The  most  likely  reason  for  such  large  differences  seems  to  be  inaccurate 
positioning  of  the  Tl  and/or  head  anatomical  coordinate  systems.  In  order 
to  analyse  this,  the  vertical  positions  of  head  and  Tl  anatomical  origins 
in  the  sled  coordinate  system  are  presented  in  Fig.  11  as  a function  of  the 
initial  neck  ler^th.  These  figures  clearly  show  that  the  neck  length  is 
strongly  related  to  the  vertical  position  of  the  Tl  anatomical  origin: 
increase  in  the  height  of  the  Tl  anatomical  coordinate  origin  causes  a 
similar  decrease  in  neck  length.  In  other  words,  the  difference  in  neck 
length  seems  to  be  mainly  due  to  errors  in  vertical  location  of  the  Tl 
anatomical  origin.  In  the  next  section  a method  will  be  proposed  to  correct 
for  this  error. 

4.4  Correction  of  Tl  coordinate  system 

The  preceding  sections  showed  errors  in  the  vertical  position  of  the  Tl 
coordinate  system  and  the  initial  Tl  orientation.  The  foilcwing  method  is 
proposed  in  order  to  correct  for  these  errors: 

- Frcm  all  tests  with  the  same  subject  a mean  value  for  the  initial  neck 
length  has  been  determined.  This  value  is  0.111  m for  subject  H0C083  and 
0.154  m for  subject  H00093. 


29 


- The  origin  of  the  new  T1  coordinate  system  is  shifted  vertically  with 
respect  to  the  laboratory  in  such  a way  that  the  initial  neck  length 
becomes  identical  to  this  mean  value. 

- The  new  T1  coordinate  system  is  initially  aligned  with  the  sled 
ratory)  coordinate  system. 

Note  that  the  orientation  of  the  T1  coordinate  system  relative  to  the  T1 
vertebra  as  defined  in  this  way,  is  dependent  on  the  impact  direction.  In 
lateral  (oblique)  direction  it  is  rotated  noninally  90  degrees  (45  degrees) 
with  respect  to  its  orientation  in  frontal  direction. 
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Fig.  10  Initial  neck  length  (i.e.  distance  between  head  and  Tl 

anatomical  origin  at  time  = 0)  in  85  tests  with  subjects 
H00083  and  H00093. 
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11  Initial  vertical  position  of  the  head  and  T1  anatauical 
origin  relative  to  the  sled  coordinate  system  as  a 
function  of  the  initial  neck  length. 
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CHAPfER  5 

ANALYSIS  OF  RELATIVE  HEAD  MOTION 


5. 1 General 

The  preceding  chapter  briefly  mentioned  the  conputer  program  which  gene- 
rates a graphical  representation  of  the  Tl  displacements  relative  to  the 
sled.  For  the  head  this  program  allows  occipital  condyles  and  anatomical 
origin  trajectories  to  be  presented  relative  to  a coordinate  systaii  with 
the  same  orientation  as  the  laboratory  (sled)  coordinate  system  and  with 
its  origin  in  the  origin  of  the  Tl  anatanical  coordinate  system.  In  addi- 
tion to  these  trajectories,  also  a trajectory  of  a point  near  the  'top'  of 
the  head  and  the  position  of  the  head  anatanical  z-axis  at  different  time 
frames  follcws  fron  this  program.  Head  twist  is  visualized  by  means  of 
small  triangles  connected  to  the  z-axis. 

Annex  B presents  for  the  30  selected  tests  the  resulting  plots.  Only  dis- 
placanents  in  the  plane  of  impact  are  presented  here  since  displacement  out 
of  this  plane  were  found  to  be  small.  For  instance  for  the  occipital  con- 
dyles and  anatanical  origin  trajectories  these  out  of  plane  displacements 
in  general  were  smaller  than  0.02  m. 

In  section  5.2  an  analysis  will  be  presented  of  these  relative  head  mo- 
tions. For  each  impact  direction,  results  of  different  tests  with  the  same 
subject  will  be  summarized  in  one  figure.  The  resulting  plots  show  that 
these  iirpact  plane  motions  can  be  approximated  quite  well  by  means  of  a 
linkage  mechanism  with  2 pivots.  Details  of  this  mechanism  will  be  dis- 
cussed in  section  5.3. 

Ihe  motion  of  this  linkage  mechanism  in  the  impact  plane  is  characterized 
by  two  degrees  of  freedan:  the  neck  link  rotation  6 and  the  head  flexion  9. 
Section  5.4  will  present  this  neck  link  rotation  as  function  of  the  head 
rotation.  In  addition  to  these  impact  plane  rotations,  in  lateral  ana 
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oblique  inpacts,  a third  significant  rotation  can  be  observed  nanely  the 
head  torsion  or  twist,  defined  as  a rotation  about  the  head  anatanical 
z-axis.  Section  5.4  will  present  this  rotation,  which  will  be  denoted  by  ij;, 
as  function  of  the  head  flexion 

5.2  Head  trajectories  and  displacements  of  the  head  anatcmical  z-axis  in 
the  impact  plane 

Fig.  12  presents  for  the  loading  phase  of  the  motion,  trajectories  in  the 
plane  of  inpact  of  the  head  anatanical  origin  and  the  occipital  condyles. 
The  position  of  the  head  anatomical  z-axis  at  three  different  time  franes 
(i.e.  50,  100  and  150  ms)  is  also  given.  The  displacements  are  presented 
relative  to  the  corrected  T1  coordinate  system  and  for  each  subject  and 
inpact  direction  the  results  are  sumnarized  in  a separate  figure.  The 
following  observations  can  be  made: 

- Trajectories  from  different  tests  are  quite  close  to  each  other  in  par- 
ticular for  the  tests  of  subject  H00093. 

- For  the  same  subject  more  severe  impact  levels  in  general  cause  larger 
head  excursions. 

- Maximum  head  excursions  in  frontal  inpacts  are  slightly  larger  than  in 
oblique  inpacts  and  much  larger  than  in  lateral  impacts. 

- In  the  initial  phase  the  head  motion  is  more  of  a translational  nature 
than  in  the  remaining  part  of  the  motion. 

- Both  the  anatanical  and  occipital  condyles  trajectories  can  be  approxi- 
mated quite  well  by  a circular  arc. 

This  last  finding  allows  the  relative  head  motions  in  the  plane  of  inpact 
to  be  approximated  in  a simple  way  by  a 2-pivot  linkage  mechanism. 
Geometrical  details  of  this  mechanism  are  provided  in  the  next  section. 

5.3  Geometrical  properties  of  the  linkage  mechanism 

The  preceding  results  suggest  that  the  relative  head  motions  in  the  inpact 
plane  can  be  approximated  quite  well  by  a 2-pivot  linkage  mechanism  (Fig. 
13).  Ihe  upper  link  of  this  mechanism  represents  the  head,  the  middle  link 
the  neck  and  the  lower  link  the  torso. 
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Fig.  12  Head  displacements  relative  to  corrected  T1  coordinate 
system  for  three  impact  directions.  Upper  curves:  head 
anatcmical  origin.  Lower  curves:  occipital  condyles. 
Stick  figures:  head  z-axis  at  50,  100  and  150  ms. 
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Fig.  13  Two-pivot  linkage  mechanism. 


The  upper  pivot  is  a joint  with  two  degrees  of  freedon  (universal  joint). 
The  first  degree  of  freedom  of  this  joint  allows  the  head  link  to  rotate 
relative  to  the  neck  link  in  the  plane  of  impact.  This  rotation  angle  will 
be  denoted  by  (J>  and  is  defined  here  as  the  angle  in  the  plane  of  impact 
between  the  z-axis  of  the  head  anatanical  and  the  corrected  T1  coordinate 
system  (Fig.  13)  (since  T1  rotations  are  neglected  in  this  study  the  z-axis 
of  the  Tl  coordinate  system  is  aligned  with  the  laboratory  z-axis).  Ttie 
mathematical  expression  for  <|)  is: 


(see  eg. (1 ) ) . 

The  second  degree  of  freedom  of  this  upper  joint  is  the  rotation  i|;  of  the 
head  about  the  head  anatomical  z-axis  indicating  the  head  torsion  or  twist. 
This  angle  will  be  derived  here  from  the  third  euler  angle  specified  by 
NBDL  (see  section  2.3).  IXie  to  definitions  of  the  coordinate  systems  the 


tan  (|)  = ^ 
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naninal  initial  value  of  the  third  euler  angle  is  0 degrees  in  frontal 
impacts,  45  degrees  in  oblique  and  90  degrees  in  lateral  impacts  (see  als*.'* 
Fig.  7).  Consequently  this  rotation  will  be  defined  by: 

= PHC03P  in  frontal  impacts 
\j;  = PHC03P-0.785  in  oblique  impacts  (7) 

^ = PHC03P-1.571  in  lateral  impacts 

The  upper  pivot  will  be  located  in  the  occipital  condylar  point.  Prelimi- 
nary results  of  a numerical  study  to  estimate  optimal  geonetrical  para- 
meters of  a 2-pivot  linkage  model  for  a set  of  6 frontal,  6 lateral  and  6 
oblique  tests  with  subjects  H00083  and  HOO093,  show  that  also  other  l(x:a- 
tions  in  the  head  can  be  chosen  for  this  pivot  (27).  However,  this  study 
did  not  shew  a significant  improvement  of  the  goodness  of  fit,  if  another 
pivot  location  was  selected. 

Ihe  Icwer  pivot  is  a pin  joint  i.e.  a joint  with  one  degree  of  freedom.  Tlie 
axis  of  the  pin  joint  is  perpendicular  to  the  plane  of  impact.  This  joint 
allows  the  neck  link  to  rotate  relative  to  the  torso  link.  Ihe  correspon- 
ding rotation  angle  will  be  denoted  by  6 and  is  defined  as  the  angle  in  the 
plane  of  impact  between  the  neck  link  and  the  z-axis  of  the  corrected  T1 
coordinate  system  (i.e.  laboratory  z-axis). 

Graphically  it  can  be  shown  that  the  circular  arc  for  approximation  of  the 
occipital  condyle  trajectories  varies  in  radius  from  0.11  m to  0.14  m.  Here 
a radius  of  0.125  m will  be  used  for  approximation  of  the  occipital  condyle 
trajectories  in  all  the  tests.  Consequently  the  neck  link  length,  i.e. 
distance  between  upper  and  lower  pivot  will  be  0.125  m.  The  corresponding 
location  of  the  lower  pivot  which  has  been  detennined  graphically  per 
subject  and  per  impact  direction  is  incorporated  in  Fig.  12.  Its  coor- 
dinates relative  to  the  corrected  Tl  coordinate  system  are  summarized  in 
Table  8. 


f 


38 


Table  8 - Location  of  lower  pivot  relative  to  corrected  T1 

coordinate  systems  for  a neck  link  length  of  0.125  m 


Impact  direction 

H00083 
X z 

(m)  (m) 

H00093 
(m)  (m) 

frontal 

-0.02 

-0.04 

-0.015 

0.00 

lateral 

-0.02 

-0.04 

-0.045 

0.00 

obi ique 

Ti.' "A  nns 

-0.025 

-0.04 

-0.035 

0.00 

Ihe  rotations  0 and  (j>  are  calculated  with  the  conputer  program  ANGLE  fron 
the  data  provided  on  the  NBDL  tapes.  Eqs.  (6)  and  (7)  show  that  (j)  and  are 
independent  of  the  geonetrical  parameters  of  the  linkage  mechanism.  Itiis  is 

not  the  case,  however,  for  the  neck  link  rotation  0.  Consequently  a change 

* 

in  pivot  location  and/or  neck  length  will  effect  the  0 value  . The  accuracy 
for  the  graphical  determination  of  the  torso  pivot  location  is  estimated  to 
be  + 0.005  m.  The  effect  of  this  on  the  accuracy  of  the  0 values  is  + 2-3 
degrees.  Numerical  techniques  are  teing  developed  that  will  allow  a more 
accurate  specification  of  link  length  and  pivot  location  (26,27). 

In  the  next  section  initial  and  maximum  values  for  the  3 rotations  0,  (j)  and 
ip  will  be  given  for  the  30  selected  tests.  In  addition  head  torsion  and 
neck  link  rotation  in  these  tests  will  be  presented  as  function  of  the  head 
flexion. 

5.4  Head  torsion  and  neck  link  rotation  as  function  c3f  head  flexion 

Values  for  the  angles  41,  4;  and  0 in  the  initial  position  (time  = 0)  will  be 
denoted  by  <|)q,  and  0Q,  respectively.  Table  9 summarizes  these  initial 
values  for  the  30  selected  tests.  The  following  observations  can  be  made: 


* In  earlier  interim  reports  and  in  ref.  (28)  and  (29)  for  seme  of 
the  tests  slightly  different  values  for  the  torso  pivot  location 
were  selected,  which  might  result  in  small  differences  in  the  0 
values  in  this  report. 
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- The  initial  head  torsion  (tj^y)  is  close  to  zero. 

- The  initial  head  flexion  in  the  plane  of  impact  (c})q)  which  is  close  to 
zero  in  lateral  impacts  shows  larger  variations  in  the  other  impact 
directions,  particularly  in  frontal  direction. 

- Ihe  initial  neck  link  rotation  in  the  plane  of  impact  (0q)  shows  a sig- 
nificant positive  value,  particularly  in  frontal  and  oblique  direction. 

The  neck  link  rotation  (0-0q)  is  presented  in  Fig.  14  as  function  of  the 
head  flexion  (4»~(|>o)*  For  each  subject  and  impact  direction  these  curves  are 
sunmarized  in  a separate  figure.  In  the  initial  phase  of  the  motion  the 
head  flexion  is  smaller  than  the  neck  link  rotation  for  all  runs  illustrat- 
ing the  presence  of  a response  which  is  frequently  perceived  as  a transla- 
tion of  the  head.  As  soon  as  the  relative  rotation  (0-0q)  - 15-30 

degrees  in  frontal  or  oblique  impacts,  changes  in  head  and  neck  angle, 
become  almost  identical.  In  other  words,  the  head  and  neck  link  are  more  or 
less  locked.  A similar  behaviour  can  be  observed  for  the  lateral  impact 
response.  However  here  the  maximum  angle  between  head  and  neck  link  is 
smaller,  namely  between  5 and  15  degrees,  while  in  addition  this  relative 
angle  appears  to  become  smaller  in  the  final  part  of  the  loading  phase. 

Table  9 summarizes  extreme  values  for  the  head  and  neck  link  rotations  in 

the  plane  of  impact,  denoted  by  d>  and  e , respectively.  Based  on  a 

^max  max 

comparison  of  (0  - ) and  (a  - in ) it  can  be  seen  that: 

max  u max  ^ 

- For  the  roost  severe  tests  in  the  three  impact  directions  head  and  neck 
link  rotations  in  frontal  flexion  are  slightly  larger  than  in  oblique 

“^impacts  and  much  larger  than  in  lateral  impacts.  Partly  this  can  be 
explained  by  the  different  impact  severities  in  the  three  impact  direc- 
tions. 

- For  similar  test  conditions  maximum  head  and  neck  link  rotations  for 
subject  H00093  are  significantly  less  than  for  subject  H00083.  This 
difference  may  be  partly  due  to  the  lower  head  mass  and  mcment  of  inertia 
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of  subject  H00093:  lioad  breadth,  head  leiKilh  and  liead  i i rcumf oi  er\co 

presented  in  table  4,  are  much  smaller  for  this  subject.  oi  o'utse  inuscde 
strength  and  the  acquired  ability  to  reduce  tlie  impact  shcx'k  al  l^'r  re- 
peated tests  may  also  contribute  to  this  difference. 


Fig.  14  Neck  link  rotation  (6  - 6^)  as  function  ol  head 
flexion  (<j>  - 4>o)  in  30  tests  with  2 siibjects. 
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Table  9 - Initial  and  maximum  values  in  degrees  for  head  flexion  head 
torsion  ^ and  neck  link  rotation  0. 


Subject 

(test 

number) 

V 

Sled 

(m/s) 

Peak 

Sled 

Acc. 

(m/s"^) 

♦o 

d> 

^max 

^ ‘^‘max 

~<l>0  ^ 

*^0  '^max 

^ 'Hrvax 

00 

A 

max 

^ ®max 
“00  ) 

FROOTAL 

H00083 

LX3530 

7.5 

49 

6 

55 

49 

9 

66 

57 

LX3536 

10.5 

73 

9 

55 

46 

- - 

- 

9 

79 

70 

LX3544 

12.4 

92 

1 

67 

66 

— — 

- 

6 

86 

80 

H00093 

LX3531 

7.6 

51 

-14 

24 

38 

- - 

- 

11 

68 

57 

LX3537 

10.1 

71 

-10 

40 

50 

- - 

- 

9 

76 

67 

LX3548 

12.2 

90 

1 

46 

45 

- - 

- 

17 

80 

63 

LX3550 

13.3 

98 

-5 

40 

45 

- - 

- 

17 

82 

65 

LX3558 

14.3 

110 

-6 

46 

52 

- - 

- 

16 

82 

66 

LX3578 

16.0 

132 

-7 

55 

62 

- - 

- 

16 

91 

75 

LX3583 

16.7 

142 

-6 

57 

63 

- - 

- 

21 

93 

72 

LX3616 

17.2 

150 

-12 

55 

67 

— — 

— 

20 

92 

72 

LATE3l^L 

HUUO'8'3"" 

LX1831 

6.5 

49 

1 

51 

50 

-2  -46 

-44 

4 

57 

53 

LX2013 

6.3 

49 

0 

40 

40 

-1  -36 

-35 

9 

51 

42 

LX2124 

3.1 

89 

-2 

35 

37 

-3  -29 

-26 

6 

39 

33 

LX2302 

6.4 

69 

-3 

50 

53 

2 -43 

-45 

3 

58 

55 

H00093 

LX2151 

3.3 

100 

0 

44 

44 

-2  -44 

-42 

3 

44 

41 

LX2326 

6.5 

60 

-3 

35 

38 

-3  -43 

-40 

3 

48 

45 

LX2355 

6.5 

50 

0 

32 

32 

-5  -39 

-34 

6 

44 

38 

OBLIQUE 

H00083 

LX2786 

7.5 

50 

5 

48 

43 

0 -14 

-14 

13 

66 

53 

LX2872 

9.0 

61 

-1 

59 

60 

3 -17 

-20 

2 

77 

75 

LX3093 

10.3 

71 

1 

40 

39 

0 -12 

-12 

6 

66 

60 

LX3102 

9.0 

60 

2 

39 

37 

1 -11 

-12 

11 

62 

51 

LX3133 

8.8 

79 

-5 

58 

63 

3 -21 

-24 

3 

83 

80 

LX3153 

8.8 

80 

5 

62 

57 

-3  -3] 

-28 

14 

93 

79 

H00093 

LX2784 

7.5 

50 

-8 

38 

38 

3 -21 

-24 

10 

68 

58 

LX2843 

3 

89 

-5 

29 

34 

-4  -25 

-21 

10 

57 

47 

LX3122 

8.9 

79 

-4 

46 

50 

-1  -29 

-28 

11 

85 

74 

LX3148 

3.9 

97 

-7 

50 

57 

-4  -31 

-27 

12 

83 

71 

LX3158 

8.8 

78 

-3 

44 

47 

-4  -29 

-25 

17 

80 

63 

LX3417 

8.8 

59 

-3 

37 

40 

-4  -26 

-22 

10 

75 

65 
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Fig,  15  presents  the  hear!  torsion  (twist)  (<}»  - ^ function  of  the 

head  flexion  ((J>  - <pQ  ) for  both  lateral  and  oblique  impacts.  No  d<ita  ior 
frontal  impacts  are  incorporated  here  since  head  tr)rsion  was  fcxmd  to  be 
negligible  in  this  impact  direction. 

For  lateral  impacts  an  almost  linear  relationship  can  be  observed  between 
head  torsion  and  head  flexion  while  the  order  of  magnitude  of  both  rota- 
tions appears  to  be  almost  identical,  particularly  for  subject  H00093. 

For  oblique  impacts  the  relation  between  head  torsion  and  head  flexion 
appears  to  be  quite  different  fron  lateral  flexion.  Peak  head  torsion 
values,  which  are  incorporated  in  Table  9,  are  significantly  less  in  obli- 
que than  in  lateral  inpacts. 
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Fig.  15  Head  torsion  ( ,{;  - ) as  function  of  liead 

((j>  - (|>q)  in  19  tests  with  2 subjects. 
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CHAPTER  6 

[joad  calculations 

6. 1 Introduction 

Itie  loads  applied  by  the  neck  to  the  head  near  the  occipital  condyles  can 
be  calculated  by  using  measurements  of  head  acceleration  and  angular  velo- 
city, if  the  head  is  regarded  as  rigid  body  and  does  not  come  into  contact 
with  any  other  object  or  body  part.  Reasons  to  perform  these  calculations 
are,  among  otherthings: 

- they  offer  an  excellent  insight  in  the  system's  behaviour  for  instance 
with  respect  to  the  role  of  the  muscle  activity 

- load-displacement  relations  can  be  used  to  formulate  dummy  performance 
requirements  and  dummy  design  specifications 

- determination  of  dynamic  properties  in  a mathematical  analog 

- it  is  generally  assumed  that  neck  loads  correlate  quite  well  with  neck 
injuries. 

The  vector  equations  for  the  force  and  the  torque  at  the  occipital  condyles 
are  obtained  by  applying  Newton's  laws: 

F = m^(ag-g)  (8) 

M = I. a + w X (I.u)  + iry^(r^^x(ag-g)  (9) 

where : 

F,  M 

I 
r 

go 

g 

a 

_g  _ 

a,  u> 


are  the  force  and  torque,  respectively  applied  by  the  neck 
to  the  head 

is  the  mass  of  the  head 

is  the  moment  of  inertia  tensor  of  the  head  (about  c.g.) 

is  the  position  of  the  head  c.g.  with  respect  to  the 

occipital  condyles 

is  the  acceleration  of  gravity 

is  the  linear  acceleration  of  the  head  c.g. 

are  the  angular  acceleration  and  angular  velocity, 
respectively,  of  the  head 
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The  linear  acceleration  at  the  head  c.g.  follows  from: 

^ ^ a X (iO) 

g a ga  ga 

where  a^  is  the  linear  acceleration  of  the  head  anatcmical  origin  and  r 

a ^ ga 

the  position  of  the  head  c.g.  with  respect  to  this  origin. 

These  equations  and  additional  relations  for  transformations  between  coor- 
dinate systems  have  been  prograimied  at  TSC  in  Cambridge,  using  as  input 
accelerations  and  velocities  derived  frcm  accelerometer  data  and  quatern- 
ians  based  on  photographic  data.  The  program  allows  expressing  the  com- 
ponents of  the  calculated  force  and  torque  with  respect  to  various  coordi- 
nate systems.  Due  to  its  three-dimensional  set-up,  this  program  can  be  used 
to  evaluate  neck  loads  in  all  impact  directions.  More  details  on  this 
program  can  be  found  in  ref.  (30). 

Results  for  the  occipital  condyle  torques  are  presented  in  section  6.2. 
This  program  was  also  used  to  estimate  the  load  at  the  center  of  the  cor- 
rected T1  origin.  Errors  made  in  these  calculations  by  neglecting  the 
inertia  effects  of  the  neck  are  expected  to  be  small. 

6.2  Occipital  condyle  and  Tl  torques 

Fig.  16,  17  and  18  present  torques  applied  to  the  head  near  the  occipital 
condyles.  Fig.  16  shows  the  ccmponent  M(j>  of  this  torque  about  an  axis 
perpendicular  to  the  plane  at  inpact  as  function  of  the  head  flexion  . 
For  each  subject  and  impact  direction  these  characteristics  are  summarized 
in  a separate  figure.  The  same  torques  are  presented  in  Fig.  17,  hcwever, 
now  as  function  of  the  relative  angle  (9"0q)  “ (<J)“<l>o  ^ between  head  and  neck 
link.  The  ccatponent  Mij;  of  the  occipital  condyle  torque  about  the  head 
anatomical  z-axis  as  function  of  the  head  torsion  (({)“(j)Q  ) is  shown  in  fig. 
18.  Results  presented  here  only  relate  to  oblique  and  lateral  impacts  since 
for  frontal  flexion  no  significant  rotations  about  the  head  anatomical 
z-axis  were  observed. 
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Frontal 
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Oblique 


Fig.  16  Torque  M near  occipital  condyles  about  an  axis 
perpendi^lar  to  the  plane  of  impact  as  function 
of  head  flexion 
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H 00083  H 000  '! ! 


Frontal 


Fig,  17  Torque  M near  occipital  condyles  about  an  axis 
perpendicular  to  the  plane  of  inpact  as  function 
of  the  relative  angle  between  head  and  neck  link 
(6-0q)  - (<|.-4.o). 
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H 00083  HOOO'n 


Fig.  18  Torque  M near  occipital  condyles  about  the  head 

anatcmicil  z-axis  as  function  of  head  torsion  (ij/njig)* 

For  frontal  impacts  the  characteristics  in  fig.  16  are  presented  along  with 
the  corridor  for  loading  proposed  by  iMertz  et  al.  (2).  The  calculated 
torque-head  flexion  characteristics  appear  to  lie  within  or  close  to  this 
corridor.  Our  calculations  seem  to  suggest  that  this  Mertz  corridor  is  too 
wide.  Largest  torque  values  (i.e.  45-55  r*n)  can  oe  observed  in  the  frontal 
tests  with  subject  H00093.  These  tests  appear  to  be  the  most  severe  tests 
in  terms  of  sled  deceleration  and  sled  velocity. 

In  lateral  impact  tests  initially  a small  positive  torque  can  be  obser- 
ved  which  is  almost  conpletely  absent  in  the  other  impact  directions.  An 
explanation  for  this  might  be  the  limited  flexibility  in  the  upper  part  of 
the  neck  in  lateral  direction  conpared  to  the  flexibility  in  the  other 
directions. 

The  occipital  condyle  torques  presented  in  fig.  17  as  function  of  the 
relative  angle  between  head  and  neck  link  confirm  the  findings  in  section 
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5.4  for  the  relative  motion  in  the  upper  pivot  of  the  linkage.  Etependent  on 
the  inpact  direction  a certain  free  range  of  motion  exists  where  the  occi- 
pital condyle  torque  is  relatively  small.  As  soon  as  this  angle  is  exceeded 
the  upper  pivot  more  or  less  gets  locked  and  the  occipital  condyle  torque 
increases. 

Fig.  18  shows  that  the  peak  value  of  the  component  M i|;  of  the  occipital 
condyle  torque  about  the  head  anatomical  z-axis  is  much  smaller  than  the 
component  of  this  torque  about  an  axis  perpendicular  to  the  irtpact  plane. 

Fig.  19  shows  an  estimation  for  the  torque  applied  to  the  head  near  the 
(corrected)  Tl  origin  as  function  of  the  neck  link  rotation  ( 0- 0q)  . Pre- 
sented is  the  component  of  this  torque  about  an  axis  perpendicular  to  the 
plane  of  inpact.  As  for  the  occipital  condyle  torques,  largest  torque 
values  appear  in  the  most  severe  frontal  impact  tests  with  subject  H00093. 
Peak  torque  values  are  lower  in  the  tests  with  subject  H00083  than  with 
H00093,  which  can  be  partly  explained  frcm  the  smaller  initial  neck  length 
in  subject  H00083  (resulting  in  a smaller  value  for  the  radius  r^^  in  eq. 
(9). 

Fig.  16-19  include  data  for  loading  as  well  as  unloading.  The  ratio  of  the 
area  between  loading  and  unloading  curve  and  the  area  between  loading  curve 
and  the  abscissa  varies  between  0.6  and  0,9  for  the  occipital  condyles  and 
Tl  torques  along  an  axis  perpendicular  to  the  impact  plane  as  function  of 
and  6,  respectively.  For  the  component  of  the  occipital  condyle  torque 
about  the  head  anatomical  z-axis  as  function  of  the  head  torsion  a much 
lower  hysteresis  is  observed. 
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Fig.  19  Torque  M near  the  corrected  T1  origin  about  an 
axis  per^ndicular  to  the  plane  of  impact  as 
function  of  neck  rotation  (e-eg)* 


■ 'va 


‘W.i 


S 


‘V. 


jL . 4|f  I ■ 15*  , 

J.M-.  .;•■  ■■'  ■ ^X' 


*c 


' > 


lEu 


—■■V. 


■H..  -jue 


r/ 


,1''^ 

fkl'  ^ 


IJr '. 

'‘•‘fcv.yb  :.h 


S«3  '•'  < 

I 15-»  , 

■ ■ ^ 


T ■-  ii. 


> ^ 


h kt«: 


U ( 


IV: 


11/^  ■;* 


'M  :<i 


K'. 


I , 


a#t»-'. 


,?>''  A;iiJ 


■■■w-  ■■W-  ■■■■;  :».?'■:■«' 


kiHt  Cm  -m.^. 

;^P  :JfV  jjllMtaa 

;»rt:^»fi- 


IP' 


:»•  m 


i!s-¥. 


/h. 


«M 


i''  , ■' ''  '']  '''’'  ‘ ’ ’ *'  fir, 

'■■  '.i  t tj 


i 


'iC-i 


4tt^ 


>iUli 


^ia> 


Tif:  ' 

». 

'*■  ',-■  v.‘  ^ 

i 


:ti 


■■■ 


V 


. -S  . 


rS 


-T?<* 


, Mm 


• I 


' h:a  :t‘t:  ’*a, 


■J 


.IS 


- ^ ■■  . 


%■'■  ■*4S^V'  !j‘'t^0>im 


r(i4j 


L'f 


/ 


•i 


^ ' fL 


[ft 


..  . ; 

‘‘‘if  -I 


-'-J 

M$:!(^)  .>W«t ' I .,«,». 


•tiUV,; 


' " rifi 

■ '.  rr"’i 


'J*! 


.<r 


GB 


1#?.. 


, w ;;i?tM,,  ,^t,,..,,  i,‘f^  , .^,, 

■^-^hfi  '-  <»»!jW>'i  'l^juiiMOJ  llbm  'te- 


s; 


. V. 


51 


CHAPTER  7 

PE3^RMANC£  REQUIREMENTS  BASED  ON  NBDL  TESTS  WITH  TWO  SUBJECTS 
7. 1 Introduction 

In  section  1.2  it  was  shown  that  performance  requirements  for  mechanical 
necks  presented  in  literature  on  this  subject  are  limited  mainly  to  the 
response  in  frontal  flexion.  Requirements  for  this  impact  direction  which 
were  developed  in  the  early  seventies,  considered  the  relation  between  the 
torque  of  the  occipital  condyles  and  the  angular  position  of  the  head. 
This,  hcwever,  is  not  a sufficient  condition  to  ensure  a humanlike  respon- 
se. 


One  method  to  define  necessary  and  sufficient  conditions  is  specification 
of  an  analog  system  with  a dynamical  behaviour  similar  to  that  of  the  real 
head-neck  system.  If  such  an  approach  is  selected,  the  analog  should  be  as 
simple  as  possible  in  order  to  limit  the  set  of  parairoters  required  to 
describe  the  analog.  The  mechanical  analog  to  be  defined  herein,  is  only 
expected  to  reproduce  the  motions  of  the  head  relative  to  the  torso  (Tl)  in 
an  impact  situation.  No  realistic  simulation  of  the  neck  behaviour  itself 
is  required.  In  fact,  since  iroasurements  have  been  conducted  only  at  the 
head  and  at  Tl,  no  actual  information  on  this  neck  behaviour  is  available. 

From  the  analysis  presented  in  chapter  5 it  follows  that  a relatively 
sinple  2-pivot  linkage  mechanism  can  be  used  to  reproduce  the  kinematic 
behaviour  of  the  head  with  respect  to  Tl.  Section  7.2  deals  with  a brief 
literature  review  of  this  type  of  analog  systems  for  the  simulation  of 
head-neck  behaviour.  In  section  7.3  joint  properties  will  be  defined  for 
the  2-pivot  linkage  on  the  basis  of  the  load  calculations  presented  in  the 
preceding  chapter.  In  this  way  together  with  mass  distribution  data,  the 
mechanical  analog  is  fully  defined. 

A second  method  to  define  performance  requirements  is  by  means  of  a set  of 
kinematic  requirements  like  displacements  and  accelerations  derived  frcra 
the  measured  test  data.  Such  a formulation  is  useful  for  instance  to  verify 
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the  response  of  a mathanatical  or  mechanical  substitute  at  a specific 
impact  level.  An  example  of  this  type  of  requirement  will  be  given  in  Vol. 
II  of  this  report  where  the  response  of  the  analog  systems  proposed  in 
section  7.3,  will  be  corpared  with  the  human  subject  response  using  mathe- 
matical simulations  . 

7. 2 Literature  review  linkage  mechanisms  for  head-neck  motion 

Analog  systems  for  the  head-neck  response  under  impact  conditions  can  be 
subdivided  into  three  categories: 

- linkage  systems  with  a limited  amount  of  degrees  of  freedom 

- mechanical  neck  system  as  used  in  various  crash  duitrny  designs 

- corplicated  mathematical  (32,33,34)  or  mechanical  (35)  models  in  which 
the  anatomical  structure  of  the  neck  (vertebra  and  ligaments)  is  repre- 
sented in  detail  (structure  models). 

This  section  deals  with  simple  linkage  systems,  while  a review  of  mecha- 
nical neck  systems  in  crash  durtny  designs  is  presented  in  ref.  (31).  Struc- 
ture models  are  too  complex  and  too  detailed  due  to  the  large  number  of 
system  parameters  involved,  to  be  used  here  for  the  specification  of  per- 
formance requirements.  This  type  of  model  is  particularly  useful  for  the 
study  of  injury  mechanisms. 

Most  linkage  systems  proposed  in  literature  for  the  head-neck  behaviour  in 
frontal  flexion,  have  two  degree  of  freedom  (e.g.  36,37,38),  although  also 
a one  degree  (39)  and  a seven  degrees  (40)  system  have  been  proposed.  The 
one  degree  of  freedom  systen,  contained  one  pin  joint  near  the  base  of  the 
neck  and  was  based  on  tests  with  human  cadavers.  In  the  seven  degrees  of 
freedom  system  seven  pin  joints  were  specified  distributed  along  the  neck 
structure.  The  pivot  locations  were  based  on  the  analysis  of  X-ray's  of 
voluntary  flexion-extension  motions  in  human  subjects  (40), 

The  two-pivot  joint  linkages  for  frontal  flexion  reported  by  Becker  (36), 
Schneider  and  EJowman  (37)  and  Frisch  and  Cooper  (38)  all  resulted  from  the 
analysis  of  NBDL  human  volunteer  tests.  The  lower  pivot  in  these  studies 
was  located  in  the  T1  origin,  which  is  different  from  the  approach  in  our 
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•Study  where  the  lower  pivot  is  located  in  the  center  o£  a circle  arc  appro- 
ximating the  occipital  condyle  trajectories.  For  the  upper  pivot  various 
locations  were  proposed  including  the  occipital  condylar  point  (38)  and  a 
point  0.04  m above  the  head  center  of  gravity  (36).  The  linkage  by 
Schtieider  and  Bowman  allows  extension  of  the  neck.  Eynamical  properties  for 
the  pin  joints  in  these  analog  systems  were  based  on  a trial  and  error 
paraaneter  fitting  process  with  a mathematical  simulation  model  until  ade- 
quate replication  of  human  response  data  was  obtained. 

For  head-neck  motions  in  lateral  impacts,  Ewing  et  al.  (20,21)  indicate 
that  this  type  of  motion  can  be  approximated  guite  well  by  a single  rota- 
tion about  a fixed  axis  in  the  laboratory.  This  axis  is  located  in  the  mid- 
sagital  plane  perpendicular  to  the  line  connecting  the  origins  of  the  T1 
and  head  anatanical  coordinate  system. 

A first  promising  attempt  to  identify  a linkage  mechanism  suitable  for 
simulation  of  frontal,  lateral  as  well  as  oblique  impacts  is  presented  by 
Becker  (41,42).  This  linkage  contains  four  degrees  of  freedon  and  is  illu- 
strated in  fig.  20.  Eight  geonetrical  parameters  specify  this  mechanism:  4 
[)arameters  that  characterize  the  position  of  head  and  Tl  with  respect  to 
the  linkage  system  and  four  parameter  that  specify  the  linkage  gecmetry 
itself.  Using  a numerical  iterative  method  these  linkage  parameters  were 
calculated  in  such  a way  that  the  linkage  approximates  photographic  data 
from  a .selected  set  of  tests  as  close  as  possible.  The  method  was  applied 
to  a .small  database  consisting  of  a voluntary  pitch  motion  and  2 lateral 
and  2 oblique  impact  tests  with  a single  volunteer  (i.e.  H00093).  Dynamical 
elements  were  not  specified  for  this  linkage. 

7 . 3 The  analog  system 

The  results  presented  in  chapter  5 indicate  that  a simple  2-pivot  linkage 
is  suitable  to  reproduce  the  observed  motions  of  the  head  relative  to  Tl. 
For  each  impact  direction  a separate  linkage  was  proposed.  It  was  shown 
that  the  upper  pivot  location  and  the  link  length  can  be  taken  identical 
for  the  3 impact  directions.  The  upper  pivot  has  one  degree  of  freedom  in 
frontal  impacts  (i.e.  head  flexion)  and  two  degrees  of  freedom  in  oblique 
and  lateral  impacts  (i.e.  head  flexion  and  head  torsion).  The  lower  pivot 
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Fig.  20  Three-dimensional  linkage  according  to  Becker  (41,42). 

has  one  degree  of  freedom  namely  the  neck  link  rotation  about  an  axis 
perpendicular  to  the  {:>lane  of  impact 

The  remaining  linkage  parameters  to  be  defined  here  are  the  lf7/^er  pivot 
hxjation,  the  initial  position  of  the  linkage,  dynamical  properties  o£  the 
pivots  and  the  head  mass  (iistribution  data.  Ihese  linkage  parameters  will 
be  determined  for  each  impact  direction  separately  on  the  basis  of  the 
results  of  the  tests  conducted  in  that  impact  direction. 
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Table  10  sunnmarizes  tlie  proposed  yeanetrical  paraneters  and  pivol,  cliaract- 
eristics.  ITae  following  additional  information  can  be  given: 
l(3wer  pivot  location:  Due  to  errors  in  the  T1  specification  it  is  not 
possible  to  specify  this  pivot  uniquely  with  respect  to  Tl  (or  torso) 
anatf^nical  landmarks.  Pivot  coordinates  relative  to  the  corrected  Tl  coor- 
dinate system,  per  subject,  were  provided  in  Table  8, 

initial  head  link  flexion  and  torsion;  Although  initial  values  for  these 
angles  are  not  completely  zero  in  the  tests,  particularly  with  respect  to 
the  initial  head  flexion,  the  deviations  from  zero  are  not  considered 
significant  to  be  accounted  for  in  the  mechanical  analog. 
initial  neck  link  rotation:  In  the  tests  significant  positive  values  for 
this  parameter  were  noted.  Ihe  values  proposed  in  Table  10  are  mean  values 
of  the  test  data. 

lower  pivot  linear  stiffness:  A linear  torque  - rotation  characteristic  is 
pro[X)sed  here  based  on  calculated  torques  near  the  Tl  origin.  For  subject 
H000B3  a correction  has  to  be  introduced  in  order  to  account  for  the  diffe- 
rence between  Tl  origin  and  lower  pivot  location.  Stiffness  data  in  Table 
10  are  based  on  a graphical  approximation  of  the  torque  - rotation  charact- 
eristics (mean  stiffness).  Ihe  accuracy  of  this  determination  is  estimated 
to  be  + 0.25  Nhi/deg.  Stiffness  in  frontal  and  oblique  direction  are  in  the 
same  order  of  magnitude,  while  in  lateral  direction  as  far  as  subject 
H00083  is  concerned,  a slightly  larger  linear  stiffness  value  can  be  ob- 
served . 

upL^er  pivot  linear  torsion  stiffness:  These  stiffness  values  appear  to  be 
relatively  small  and  are,  just  as  for  the  lower  pivot  stiffness  based  on  a 
graphical  estimation  (accuracy  = 0.1  Nm/deg.). 

upper  pivot  free  range  of  motion  in  the  plane  of  impact;  This  parameter  can 
be  determined  graphically  fron  the  torque-rotation  characteristics  provided 
in  fig.  17.  For  frontal  and  oblique  impacts  the  range  of  motion  is  of  the 
same  of  order  of  magnitude;  for  impacts  in  lateral  direction  a smaller 
value  is  found.  Accuracy  for  this  quantity  is  + 3 degrees.  No  stiffness 
value  for  the  joint  stop  (locking  mechanism)  can  be  specified  here  on  the 
basis  of  the  available  data. 

In  addition  to  the  preceding  parameters  the  mass  distribution  of  the  head 
has  to  be  defined.  In  this  study  this  will  be  based  on  the  mass  distribu- 
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tiun  of  subject  H00083  and  H00093  by  defining  average  values  from  the  mass 
distribution  data  provided  in  table  6.  For  the  neck  link  the  mass  will  be 
taken  zero  since  the  possible  inertia  effects  of  the  neck  were  neglected  in 
the  T1  torque  calculations. 


Table  10  - Linkage  parameters  for  the  loading  phase  based  on  30 
impact  tests  with  two  subjects  (H00083  and  H00093) 


Parameter 

GP^CrRICA^[.  PARAMETERS 
upper  pivot  location 

link  length 
lower  pivot  location 
initial  head  link  flex, 
initial  head  link  tors, 
initial  neck  link  rot. 


Frontal 

direction 


occipita  L 
condyles 
0.125  m 
see  table  8 
0 
0 

14  deg. 


DYNAMICAL  PARAMETERS 
lower  pivot  linear 

stiffnes  1.25  Nm/deg. 

upper  pivot  linear 

torsion  stiffness 

upper  pivot  free  range 

of  motion  in  the  plane  ^ 

of  impact  25 


Lateral  Oblique 

direction  direction 


occipital 
condyles 
0. 125  m 
see  table  8 
0 
0 

5 deg. 


occipital 
condyles 
0.125  m 
see  table  8 
0 
0 

10  deg. 


1 . 5 Nm/deg . 1.25  Nm/deg . 

0.25  Nm/deg.  0.5  Nm/deg. 
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CHAPTER  8 

ANALYSIS  OF  ADDITIONAL  TESTS  WITH  OrHEK  SUBJfciCrS 


8. 1 Introduction 

In  addition  to  the  tests  with  subject  H00083  and  H00093  discussed  before, 
results  will  be  presented  here  of  a number  of  earlier  frontal  and  lateral 
tests  with  other  human  subjects.  The  frontal  tests  were  conducted  in  the 
period  between  1967  and  1969  at  Wayne  State  University  and  are  reported  by 
NBDL  in  reference  (7,16).  No  experimental  results  on  magnetic  tape  were 
available  for  these  tests,  therefore  only  a rough  analysis  of  these  tests 
could  be  made  (see  section  8.2). 

The  lateral  tests  to  be  considered  here  were  reported  by  NBDL  in  reference 
(10).  In  contrast  to  the  frontal  tests,  for  this  lateral  test  series,  data 
tapes  with  experimental  results  were  available.  An  analysis  of  these  tests 
was  presented  at  the  27th  Stapp  Conference  (28)  together  with  the  lateral 
impact  tests  with  subject  H00083  and  H00093.  Section  8.3  gives  a brief 
sunmary  of  the  findings  of  this  analysis. 

8. 2 Frontal  impact  tests 

A detailed  description  including  test  method,  subject  anthropometry  and 
test  results  of  the  additional  frontal  tests  to  be  considered  here,  is 
given  by  Ewing  and  Thonas  (7,16).  These  tests  were  conducted  at  Wayne  State 
University.  Test  conditions  were  slightly  different  fron  those  of  the  tests 
with  subjects  H00083  and  H00093.  For  instance,  the  sled  acceleration  - time 
history  was  triangular  rather  than  trapezoidal  (see  Fig.  21).  Results  of 
the  most  severe  tests  with  10  subjects  are  used  here.  Table  11  summarizes 
the  most  significant  test  conditions. 

Table  12  summarizes  the  type  of  information  available  in  reference  (7) 
which  was  used  for  our  analysis.  From  this  data  it  was  possible  to  recon- 
struct the  location  of  the  occipital  condyles  at  certain  points  in  time 
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relative  to  the  sled  coordinate  system.  The  shape  of  the  center  of  gravity 
trajectories  was  found  to  be  circular  and  fron  the  graphical  reconstructed 
occipital  condyles  it  followed  that  for  most  of  the  subjects  a 2 pivot 
linkage  system  with  a link  length  of  0.125  m can  be  used  to  approximate  the 
observed  motions. 


Fig.  21  Conparison  of  sled  accelerations. 


Table  11  - Test  characteristics  of  earlier  tests  with 

10  different  subjects  in  frontal  impacts  (7). 


Subject 

Sled 

Peak 

Rate  of 

velocity 

sled 

onset 

change 

3CC  • 2 

3 

(m/s) 

(m/s  ) 

(m/s  ) 

003 

13 

95 

3500 

004 

13 

94 

3300 

007 

14 

96 

6300 

009 

13 

90 

3300 

010 

13 

89 

3300 

Oil 

12.5 

89 

3800 

012 

13.5 

98 

5300 

013 

12.5 

91 

5300 

015 

14 

99 

4600 

016 

12 

84 

3400 

Table  13  summarizes  the  initial  and  extreme  values  for  the  head  flexion 
which  could  be  determined  simply  from  the  information  available  for  the 
angular  displacements  (PESi).  Initial  and  maximum  values  for  the  link  ro- 
tation were  graphically  estimated  from  the  reconstructed  occipital  condyle 
locations. 
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r.ibli>  \2  - Scb'cl  ion  oi  t:c\sL  tosu  1 1 s provided  in  ret;er(?noe  (7). 


oC  IJahn 

Center  ot  gravity  trajectory 


Tl  angular  displacement 


Head  angular  displacement 


Center  of  gravity  accele- 
ration 


fiead  angular  acceleration 


[Adscript  ion  

Coordinate  of  center  of  gravity 
(MCX  and  MCZ)  relative  to  Tl 
coordinate  system 

Angular  displacement  (NRN)  of  Tl 
coordinate  system  as  function  of 
time  relative  to  laboratory 
coordinate  system 

Angular  displacement  (PDH  and  PHN, 
respectively)  as  function  of  time 
relative  to  laboratory  and  Tl 
coordinate  system,  respectively. 

Linear  accelerations  (4XM  and 
4ZM)  as  function  of  time  relative 
to  laboratory  x and  z axis 

Angular  acceleration  (RAH)  as 
function  of  time 


Table  13  - Initial  and  maximum  values  for  the  head  rotation  ^ and 
the  link  rotation  0 (in  degrees). 


Subject 

max 

(4'  - 4>o) 

max  “ 

®0 

6 

max 

(6  -Bq) 

max 

003 

3 

81 

78 

36 

90 

54 

004 

- 7 

80 

87 

34 

98 

64 

007 

- 3 

90 

93 

10 

no 

100 

009 

- 6 

71 

77 

27 

102 

75 

010 

-15 

71 

86 

19 

105 

86 

Oil 

0 

72 

72 

32 

108 

76 

012 

-13 

72 

85 

9 

112 

103 

013 

0 

76 

76 

28 

LOO 

72 

015 

3 

80 

77 

28 

104 

76 

016 

3 

71 

68 

34 

103 

69 

■ITie  accuracy  for  these  graphically  determined  neck  link  rotations  is  esti- 
mated to  be  + 3 degrees.  Just  as  for  the  tests  with  subject  H00083  and 
H00093,  a significant  positive  value  for  the  initial  neck  link  rotation  can 

be  observed.  Maximum  neck  link  rotations  (9  ~9„)  are  of  the  same  order  of 

max  0 

magnitude  as  for  the  test  with  subjects  H00083  and  H00093,  while  values  for 

the  head  flexion  (a  -An)  appear  to  be  larger  for  most  of  these  tests. 

anax 


Peak  torque  values  near  the  occipital  condyles  for  the  1973  datalvise  were 
reported  by  Ewing  anil  Thomas  (16)  to  vary  from  30  Nm  to  50  Nm,  which  is 
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close  to  the  peak  values  in  Fig.  16.  In  contrast  to  our  findings,  Swing  and 
ThOTas  data  showed  significant  pcxsitive  torque  values  in  the  initial  phase 
of  the  motion.  This  difference  can  be  partly  explained  by  different  values 
used  by  Ewing  and  Thomas  for  the  center  of  gravity  location.  Our  calcula- 
tions were  based  on  more  recent  cadaver  data  (25)  which  showed  a larger 
distance  of  the  head  c.g.  to  the  anatomical  origin  than  in  the  data  used  by 
P]wing  and  Thomas.  Moreover,  larger  initial  angular  accelerations  were 
measured  in  the  Ewing  and  Thcmas  tests.  The  reasons  for  these  larger  values 
are  not  fully  clear  yet,  but  the  fact  that  in  these  earlier  tests  slightly 
different  test  conditions  were  used  might  contribute  to  these  differences. 
A more  detailed  discussion  on  these  aspects  is  given  by  Spenny  and  Wismans 
(43). 

Fran  the  mass  distribution  data  provided  in  ref.  (16)  and  the  linear  and 
anjular  accelerations  given  in  ref.  (7)  a rough  estimate  could  be  made  for 
the  maximum  torque  near  the  center  of  rotation  in  the  torso.  These  torques 
were  found  to  vary  from  59  Nra  to  114  Nm  which  is  close  to  the  peak  values 
in  fig  19. 

8.3  Lateral  impact  tests 

The  lateral  tests  which  will  be  considered  here  were  conducted  by  NBDL  for 
NHTSA  under  contract  DOT-HS-4-0852.  In  total  34  runs  with  6 human  volun- 
teers and  12  runs  with  2 chimpanzees  were  conducted.  Photographic  data  on 
tape  for  the  chimpanzee  tests  were  limited  to  a time  frame  before  onset  of 
the  sled  acceleration,  so  that  no  further  analysis  of  these  animal  tests 
was  made.  Out  of  the  34  tests  with  6 human  subjects  a subset  was  defined  of 
9 of  the  most  severe  tests  with  four  different  subjects.  Table  14  summa- 
rizes these  tests  and  the  most  important  test  conditions.  The  impact  seve- 
rity in  this  database  is  of  the  same  order  of  magnitude  as  the  impact 
severity  in  the  H00083  and  H00093  lateral  flexion  tests. 


61 


Table  14  - Test  Characteristics  of  16  selected  runs  witli  6 human 
volunteers  in  lateral  flexion 


Subject 

Run 

Sled 

Peak 

Rate  of 

Number 

Velocity 

Sled 

Onset* 

Change 

Acc. « 

(m/s) 

(m/s^) 

(m/s"^) 

H00044 

LX1504 

6.6 

60 

1246 

H00U44 

1^(1512 

6.2 

69 

1519 

H00044 

LX1528 

6.3 

74 

1644 

H00064 

LX1507 

6.5 

59 

1202 

H00064 

LX1513 

6.2 

69 

1506 

H00065 

LX1505 

6.6 

61 

1233 

H00065 

LX1510 

6.  2 

69 

1526 

400067 

LX1509 

6.2 

68 

1498 

400067 

LX1525 

6.2 

74 

1708 

* Slope  of 

the  best 

least 

square  line 

fit  of 

the  rising  portion 

of  the  sled  acceleration  profile 

between  20%  and 

50%  of  the 

peak  sled 

1 acceleration. 

Table  15  - 

Initial  and  maximum  head 

flexion  <j) 

and  1 

neck  link 

rotation 

0 (in  degrees) 

in 

lateral 

flexion 

Subject  nr. 
(test  nr. ) 

♦o 

<t> 

max 

'I* 

max  ^ 

0 

max 

0 -0Q 

max  “ 

H00044 

1504 

-1 

39 

40 

0 

47 

47 

1512 

6 

53 

47 

5 

57 

52 

1528 

3 

48 

45 

5 

56 

51 

H00064 

1507 

3 

61 

58 

-4 

61 

65 

1513 

4 

64 

60 

2 

65 

63 

400065 

1505 

1 

59 

58 

7 

67 

60 

1510 

6 

63 

60 

7 

67 

60 

H00067 

1509 

4 

59 

55 

-1 

58 

59 

1525 

6 

60 

64 

4 

62 

58 

Results  of  these  tests  were  found  to  agree  f^ite  well  with  the  test  results 
obtained  with  subjects  H00083  arid  H00093.  A similar  linkage  as  proposed  in 
chapter  7 for  the  lateral  impacts  could  be  used  to  characterize  the  obser- 
ved displactsments,  namely  a 2-pivot  1 inkage  mechanism  with  the  upper  pivot 
located  at  the  occipital  condyles  and  a link  length  of  0.125  m.  The  lower 
pivot  location  in  the  torso  as  for  subject  H00083  and  H00093  was  found  to 
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Uj  diCferent  for  each  subject.  A rough  estimation  tor  these  lotraLious 
resulted  in  a distance  varying  fnxn  0 - 0.02  m to  the  mici-sagittal  plane, 
in  the  direction  of  the  thrust  vector  of  the  accelerator. 

In  order  to  illustrate  the  similarity  in  kinematics  a comparison  is  pre- 
sented in  Fig.  22  for  the  head  torsion  (i|;-i|»o)  neck  link  rotation  {6-8o) 
as  function  of  head  flexion  (<j>-(j))*.  A summary  of  extreme  value  for  the  head 
flexion  <j,  and  the  neck  link  rotation  e is  given  in  table  15.  Maximiim  rota- 
tions ((j)~(j)Q)  and  (0-9q)  appear  to  be  slightly  larger  for  some  of  the  sub- 
jects in  this  database  compared  to  the  tests  with  subject  H00083  and 
H00093. 

For  the  loads  at  the  occipital  condyles  and  at  the  torso  center  of  rotation 
the  tests  in  both  databases  showed  good  agreement  (28). 


* e values  in  this  figure  for  subject  H00083  and  H00093  may  slightly  vary 
from  data  presented  in  chapter  5,  due  to  the  use  of  a different  method 
used  for  Tl  coordinate  system  correction  and  for  specification  of  the 
lower  pivot  location.  The  torso  pivot  was  taken  identical  for  all  6 
subjects  based  on  an  average  value  for  all  6 subjects  (i.e.  0.02  m fran 
the  mid-sagittal  plane). 


r 
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-(Pg  (rad) 


Fig.  22  Neck  link  flexion  {0-0q)  and  head  link  torsion  ('I'-'J'q) 
as  function  of  head  flexion  («|>-(j)Q)  in  16  lateral 
impacts  with  subjects  H00044,  fl00064,  H00065,  H00067, 
H000B3  and  H00093. 
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Cl-lAPrER  9 

DISCUSSION  AND  CONCUJSIONS 


'Ihe  study  presented  here  is  part  of  a research  program  to  develop  a mecha- 
nical neck  with  omnidirectional  biofidelity.  An  analysis  is  presented  of  30 
tests  with  2 human  volunteers  exposed  to  frontal,  oblique  and  lateral 
impacts. 

The  most  important  finding  of  this  study  is  that  the  observed  human  head- 
neck  response  can  be  represented  adequately  by  a linkage  system  with  2- 
pivots.  For  the  three  impact  directions  a separate  mechanical  analotj  is 
proposed.  For  frontal  impact  this  system  has  2 and  for  lateral  and  oblique 
impacts  3 degrees  of  freedom.  These  degrees  of  freedom  are  a neck  link 
rotation  in  the  plane  of  impact,  a head  rotation  in  the  plane  of  impact  and 
a head  torsion  about  the  head  anatomical  z-axis  respectively.  This  last 
de<jrfie  of  freedom  is  absent  in  frontal  impacts.  The  neck  link  length  and 
the  upper  pivot  location  roughly  can  be  taken  identical  in  all  three  impact 
directions.  For  sane  of  the  ranaining  parameters  of  the  mechanical  analog 
differences  for  the  three  impact  directions  can  be  observed  (see  7.3).  In 
Vol.  II  of  this  report  the  proposed  mechanical  analog  systems  will  be 
simulated  mathematically  and  the  model  response  will  be  compared  with  the 
human  volunteer  behaviour. 

In  this  study  separate  systems  are  proposed  for  the  three  impact  direc- 
tions. As  an  alternative,  one  single  three  dimensional  linkage  could  also 
have  been  proposed  to  simulate  all  three  inpact  directions.  Such  a mecha- 
nis-n  would  p<r5sess  5 degrees  of  freedom:  2 degrees  of  freedom  in  the  lower 
pivot  and  three  degrees  of  freedom  in  the  upper  pivot.  In  such  a system  the 
neck  link  rotation  is  not  limited  to  a planar  motion  in  the  plane  of  im- 
pact. The  major  problem  in  defining  such  a system  is  the  initial  orienta- 
tion of  the  neck  link  in  lateral  flexion.  If  this  initial  orientation  (i.e. 
5 degrees)  can  be  neglected,  fonnuLatlon  of  such  a 3-dimensional  system  is 
not  expected  to  cause  significant  problems.  Attractive  aspects  of  such  a 
three-dimensional  system  are,  among  other  things,  the  possibility  to  pre 
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diet  the  head  response  in  other  (Jirections  than  frontal,  lateral  nc  oblique 
directions  and  the  possibility  to  d<itermine  out  of  impact  plane  motions 
which  were  neglected  in  the  [present  analysis. 

Q-jornetrical  linkage  parameters  in  this  study  have  been  estimated  using 
simple  graphical  techniques.  In  addition,  numerical  techniques  were  deve- 
loped and  applied  to  a limited  set  of  tests  which  will  allow  a more  accu- 
rate specification  of  these  parameters  (27).  Kspecially  for  a database  with 
a large  number  of  human  volunteer  and  caciaver  tests  a numerical  approadi  is 
expected  to  be  beneficial. 

Tne  head  link  rotation  (flexion  and  torsion)  can  be  detennined  directly 
from  the  test  data  provided  by  NBDL.  However,  this  is  tv^t  the  case  for  the 
neck  link  rotation.  A different  correction  method  for  the  Tl  coordinate 
system  for  instance  and  the  limited  accuracy  of  the  graphical  pivot  speci- 
fication will  effect  the  0 values. 

One  of  the  most  significant  parameters  in  the  mechanical  analog  is  the 
lower  pivot  stiffness.  This  stiffness  was  estimated  from  torque  data  near 
Tl,  rather  than  based  directly  on  torque  calculations  at  the  lower  pivot 
location.  It  is  strongly  recoranended  in  subsequent  test  series  to  calculate 
this  torque  at  the  location  of  the  lower  pivot  in  order  to  improve  the 
accuracy  of  the  lower  pivot  stiffness  specification  (see  also  Vol.  IT). 

Tn  this  study  pivot  stiffness  data  are  derived  graphically  from  torque- 
rotation  characteristics.  The  limited  accuracy  of  this  graphical  approach 
should  improved  by  means  of  numerical  (least  squares)  approximation 
techniques. 

Ihe  findings  in  this  study  are  primarily  based  on  30  tests  with  two  sub- 
jects. These  findings  were  roughly  compared  with  some  of  the  earlier  tests 
by  NBDL  (see  chapter  8)  and  in  general  good  agreement  could  be  observed.  It 
is  strongly  reconrmended  that  additional  tests  with  other  human  subjects  are 
analysed  to  verify  in  more  detail  the  findings  presented  here. 
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Hie  NBDL  tests  were  conducted  wi.tfi  u variety  of  test  conditions  for  rela- 
tively lf>w  j-levels.  The  severity  of  the  most  severe  human  volunteer  test 
(i.e.  15.3  g,  17,2  m/s)  is  still  less  severe  than  the  standard  Part  572 
neck  pendulum  calibration  test,  while  on  the  other  hand  this  calibration 
test  is  much  less  severe  than  for  instance  a 35  mph  car  crash  test  envi- 
ronment. Since  the  pro[>ised  analog  is  only  valid  for  low  severity  impacts, 
i.e.  the  NBDL  human  volunteer  test  conditions,  additional  information 
should  be  obtained  for  higher  exposure  levels  frcan  human  cadaver  tests.  If 
such  data  liecome  available,  adjustment  of  the  proposed  analog  systems  might 
be  necessary.  In  conjunction  with  such  tests  mathematical  simulation  could 
be  conducted  in  order  to  correct  for  the  absence  of  muscle  activity. 

The  mechanical  analog  systems,  as  defined  herein,  implicitly  constitute  a 
[jerformance  requirement.  It  is  also  possible  to  explicitly  specify  the 
(v-cessary  and  sufficient  conditions  for  achieving  the  major  response  by 
miMns  of  kinenatic  requirements  (displacements,  accelerations)  for  several 
iiTpact  levels.  Such  a formulation  is  desirable  for  establishing  a standard 
t.o  b»e  used  in  the  laboratory  for  verifying  response  of  a duniny.  In  such  a 
formulation  the  measured  volunteer  test  data  is  the  performance  require- 
ment, thereby  eliminating  the  loss  of  rigor  due  to  rigid  body  modelling 
assumptions  and  due  to  interpretation  in  formulating  the  mechanical  sub- 
stitute's characteristics  from  the  test  data.  An  example  of  this  type  .of 
requirements  is  provided  in  Vol  II  of  this  final  report  where  the  response 
of  an  analog  system  which  is  simulated  mathematically,  will  be  compared 
directly  with  the  human  volunteer  response. 
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APPENDIX  A 

(XMPARISON  OF  SLED  ACCELERATION  AND  HORIi^ONTAL  'IT 
ACCELERATION  IN  30  SELECTED  TESTS  WITH  SUBJECT 
H00083  AND  H00093, 
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Annnex  A 


Comparison  of  sled  acceleration  and  horizontal  T1  acceleration  in 
30  selected  tests  with  subject  H00083  and  H00093. 
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APPENDIX  B 


RELATIVE  HEAD  MOTIONS  IN 
30  SELECTED  TESTS  IN  PLANE  OF  IMPACT 
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